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ABSTRACT 

Using semi-analytic models of galaxy formation set within the Cold Dark Matter 
(CDM) merging hierarchy, we investigate several scenarios for the nature of the high- 
redshift (z ;> 2) Lyman-break galaxies (LBGs). We consider a "collisional starburst" 
model in which bursts of star formation are triggered by galaxy-galaxy mergers, and 
find that a significant fraction of LBGs are predicted to be starbursts. This model 
reproduces the observed comoving number density of bright LBGs as a function of 
redshift and the observed luminosity function at z ~ 3 and z ~ 4, with a reasonable 
amount of dust extinction. Model galaxies at z ~ 3 have star formation rates, half-light 
radii, I— K colours, and internal velocity dispersions that are in good agreement with 
the data. Global quantities such as the star formation rate density and cold gas and 
metal content of the Universe as a function of redshift also agree well. Two "quiescent" 
models without starbursts are also investigated. In one, the star formation efficiency in 
galaxies remains constant with redshift, while in the other, it scales inversely with disc 
dynamical time, and thus increases rapidly with redshift. The first quiescent model is 
strongly ruled out as it does not produce enough high redshift galaxies once realistic 
dust extinction is accounted for. The second quiescent model fits marginally, but un- 
derproduces cold gas and very bright galaxies at high redshift. A general conclusion 
is that star formation at high redshift must be more efficient than locally. The colli- 
sional starburst model appears to accomplish this naturally without violating other 
observational constraints. 

Key words: galaxies: high redshift - galaxies: formation - galaxies: evolution - 
galaxies: starburst - cosmology: theory 



1 INTRODUCTION 

With the dramatic recent advances in observational astron- 
omy, more and more pieces of the puzzle of galaxy formation 
and evolution are becoming available. Some of the important 
puzzle pieces include the number densities, colours, sizes, 



niques were exploited by Madau et al. (1996 ) to identify 
high-redshift candidate s in the Hubble Deep Field (HDF) 
(Williams et al. 1996). Extensive spectroscopic follow-up 
work at the Keck telescope has veri fied the accuracy of 
the photometric selection technique (Stcidel et al. 1996b; 



morphc logies, internal velocity dispersions and star forma- 



Lowenthal et al. 1997; Adelberger et al. 199£; Steidel et al 



tion ra tes of bright star-forming galaxies spanning a redshift 
range Lorn z ~ to z ~ 6, and the complementary infor- 



199£ ). The morphologies and siz es of these objects have been 
studied using the HDF sample ([Giavalisco, Steidel, fc Mac" " 



chetto 1996 



Lowenthal et al. 1997), and their clustering 



mation on the neutral hydrogen and metal content of the properties have been measured using the growing sample 
Universe to z ~ 4 obtained from quasar absorption systems. of hundreds (now approaching a thousand) of Lyman-break 
However, it still remains to fit these pieces together into a galaxies (LBGs) w i th spectroscopically confirmed redshifts 
compre hensive and compelling theoretical framework. (feteidel et al. 199£ ; Giavalisco et al. 1998| ; Adelberger et al 



OiL winHnw nntn thp high rpHghift (y ^> 9) TTniwrap Viag 

been expanded tremendously by the "Lyman-break" photo- 
metric s election technique developed by Steidel and ggU ab 



orators 



199S ) . Simila r techniques have b een used to identify galax- 
ies at z ~ 4 (Steidel et al. 199£ ) , and a handful of o bjects 
have been discov ered with confirmed redshifts z J> 5 (Stern 



fc Spinrad 1999) 



Steidel et al. 1996c ). This technique uses specially developed 
filters to exploit the redshifted Lyman-limit discontinuity 
in order to identify high-redshift candidates. Similar tech- 



Soon after the discovery of Lyman-break galaxies, var- 
ious inte rpretations of them w ere proposed. A vi ew put for- 
ward by Steidel et al. (1996a ) and reiterated by |Giavalisco" 



© 0000 



RA£ 



2 R. S. Somerville, J. R. Primack & S. M. Faber 



Steidel, fc Macchet to (1996), |Steidel et al. (199^ ), and |Adel- 
berger ;t al. (1998), is that LBGs are located in the centres 



of massive dark matter haloes (M :> 10 Mq) and have 
been forming stars at moderate rates over a fairly long time 
scale (> 1 Gyr). This scenario supposes that LBGs are the 
direct progenitors of today's massive, luminous ellipticals 
and spheroids and that such objects had almost completely 
assembled and begun f orming stars at modera te and fairly 
constant rates by z ~ 3. Lowenthal et al. (1997) also equated 



limit by short, intense bursts of star formation. Trager et aL 
(1997|) further proposed that these small objects would later 
be tidally destroyed as they fell into the potential well of 
larger galaxies to produce a Population II stellar halo like 
that in the Milky Way. 

The impressive body of observations at high redshift 
would seem to offer a tantalizing promise to provide impor- 
tant constraints on theories of galaxy formation. As usual, 
the difficulty lies in making a connection between the vis- 
ible but ill-understood portions of galaxies and the dark 
matter haloes that are hidden from view but well modelled 
within the Cold Dark Matter (CDM) framework. Building 
on the above massive/early-assembly picture for LBGs, sev- 
eral workers investigated a model in which each LBG is 



simply changing this ingredient can dramatically change our 
predictions about the high redshift Universe, and discuss 
which — if any — of the recipes can be ruled out by com- 
parison with the observational data. 

We will define "quiescent" star formation to be the stan- 
dard mode of star formation that occurs within galactic discs 
whenever cold gas is present. It is the dominant mode of star 
formation in the present-day Universe, and there is observa- 
tional evidence that its efficiency in nearby galaxies is related 
LBGs with today's massive spheroids but suggested instead to the internal properties of galactic discs, appare ntly either 

that L BGs are still actively assembling and that many ex- to th e surface density or to the dynamical time ( Kennicutt 

ist in s [nall-mass haloes raised briefly above the detection 199£ ). Since both of these properties are slowly varying, the 

star formation rate under quiescent star formation changes 
only slowly with time. 

We consider also a second mode of star formation in 
which stars are created with dramatically increased effi- 
ciency over relatively short time scales, termed "starbursts." 
Various physical mechanisms might produce such bursts, 
but here we will assume that they are triggered by galaxy- 
galaxy mergers. We term these "collisional" starbursts to 
distinguish them from bursts triggered by other means. Sub- 
stantial observational and theoretical evidence exists in sup- 
port of the collisional starburst phenomenon. For example, a 
strong correlation is observed b etween starburst activity and 
interactions in local galaxies (Kennicutt 1998), and high- 



associated with a dark matter halo and luminosit^scales 



tightly 



resolution N-body simulations including gas dynamics ( M 
hos fc Hcrnquist 1994 Miho-s fc Hernquist 1995| ; |Mihos"~& 



et al. l piil ; [jing fc Suto 1998| ; |Wechsler et al. 1998| ). This 



with lia.krTmsTT^ ^^^Tugrh^lT)!)!!^ Hcrnquist 1996 Barnes fc Hcrnquist 19!)(i| ) have shown that 



work has demonstrated that, in such a picture, the observed 
number density and clustering of LBGs at z ~ 3 can be 
simultaneously accounted for if there is approximately one 
galaxy per massive halo (M > 10 12 M@, where this mass de- 
pends somewhat on the cosmological model) . We shall refer 
to this class of simple models as "massive halo" models. It 
arises quite generically from such models that LBGs must be 
substantially more clustered than the underlying dark mat- 
ter — in short, bright galaxies at high redshift are highly 
biased. This prediction has been em phasized as a major suc - 
cess of the massive halo picture (e.g., Adelberger et al. 1998). 

However, the sharp lower mass cutoff and tight link be- 
tween halo mass and luminosity in massive halo models are 
doubtless too simplistic. A step forward in realism is pro- 
vided by semi-analytic models, in which the infall rate of 
cooling gas determines the star formation rate, and hence 
the luminosity. Both are modelled as a function of halo 
mass, which grows via hierarchical clustering, accompanied 
by a sch ematic law for gas c ooling as a function of halo gas 
density. Baugh et al. (19*9^ , hereafter BCFL) showed that 
semi-analytic models are broadly consistent with the simple 
massive halo picture, and that they also match the obs erved 
abund ance and cl ustering of LBGs at z ~ 3 (see also Gov- 
ernato fct al. 1998[). 



While these results are encouraging, there is an impor- 
tant caveat. Agreement with the observed number densities 



collisions can trigger a bar-like instability that efficiently 
drives gas into centrally concentrated, high density knots, 
creating the conditions that are likely to result in a star- 
burst. 

Our terminology implies a dichotomy between "quies- 
cent" and "bursting" star forma tion that may be somewhat 
artificial — the recent work of Kennicutt (199S ) suggests 
that star formation in both quiescent (normal) and star- 
burst galaxies has the same dependence on the local gas 
density. By funneling gas into a small central region, the 
interaction may simply produce the elevated gas densities 
that in turn lead to enhanced levels of star formation. Still, 
in interpreting the high redshift observations it is impor- 
tant conceptually to determine what kind of process actu- 
ally dominates early star formation; the difference between 
internally-governed versus externally-triggered star forma- 
tion is highly significant. For example, it determines whether 
we interpret the number and luminosity density of high red- 
shift galaxies as reflecting basically their masses and internal 
properties, or instead the merger rate at that epoch and the 
efficiency of the gas inflows produced by these mergers. 

In this paper, we investigate a scenario in which most 
of the observable LBGs are collisional starbursts. Despite 
the fact that such collisio ns inevitably arise in hierarchi- 
cal clustering models (e.g., Kolatt et al. 2000), small-object 
collisions would not have been resolv ed in existing hydro- 
dynamical simulations (for example Katz, Hernquist, & 



of LBG i is achieved in the above models only if internal dust Weinberg 1999), while previous semi-analytic work may 



extinction in LBGs is ignored. Much evidence has now accu- 
mulated to the effec t th at dust in LBGs is in fact highly sig- 
nificant (see Section 2.4). The shortfall in LBG numbers that 
results when realistic dust extinction is included (which we 
will quantify) motivates us to explore other approaches, in 
particular, other modes of star formation. We will show how 



have also underestimated the importance of collisional star- 
bursts because of their assumed physical recipes. We inves- 
tigate whether the numbers and properties of high-redshift 
starbursts are compatible with observations of LBGs, and 
whether this scenario leads to self-consistent agreement with 
global quantities such as the star formation rate density and 
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the cold gas and metal content of the Universe as a function 
of redshift. We also consider models containing only "quies- 
cent" star formation and discuss whether the collisional star- 
burst picture is merely compatible with the data or whether 
some contribution from bursts seems to be required. Our 
conclusion is that quiescent models can marginally match 
all data provided that early star formation is highly accel- 
erated, but that burst models fit all data naturally and are 
therefore preferred. Either way, we conclude that early star 
formation must be much more efficient than locally. 

The outline of this paper is as follows. In Section ^, we 
give a brief introduction to the semi-analytic models, includ- 
ing our simple approach for including collisional starbursts. 
In Section tl we present the predictions of our fiducial mod- 
els and compare them with the observations. In Section W, 
we compare our results with previous work, and in Section g 
we discuss the sensitivity of our results to various assump- 
tions, including the cosmology, stellar population synthesis 
and dust modelling. We summarize our results and conclude 
in Section [| The non-specialist reader may wish to skip Sec- 
tions ^ and which are somewhat technical. The main re- 
sults of these sections are summarized in a general way in 
Section |^. 



2 SEMI-ANALYTIC MODELLING 
2.1 Basics 

We use semi-analytic techniques to model the formation and 
evolution of galaxies in a hierarchical clustering framework. 
Our models include the effects of gravitation on the for- 
mation and merging of dark matter haloes, the hydrody- 
namics of cooling, star formation, supernovae feedback and 
metal production, galaxy-galaxy merging, and the evolution 
of stellar populations. Our Monte Carlo-based appr oach is 



similar in spirit to the models origin ally presented by Kauff- 



mann , White, fc Guiderdoni (1993 , hereafter KWG93) and 



Cole et al. (1994 hereafter CAFNZ94), and subsequently de- 
veloped by these groups (hereafter referred to as the "Mu- 
nich" and "Durham" groups) and others in numerous pa- 
pers. The semi-a naly tic models used here are de scribed in 
Somerville (1997) and Somerville & Primack (1999, hereafter 
SP). As shown in SP, these models produce good agreement 
with a broad range of local galaxy observations, including 
the Tully-Fisher relation, the B and K-band luminosity func- 
tions, cold gas contents, metallicities, and colours. We refer 
the reader to SP for a more comprehensive review of the 
literature and a more detailed description of our models. 
Below we sketch our approach briefly. 

Th e framework of the semi-analytic approach is the 



merging history" of a dark matter halo of a given mass, 



identified at z = or any other redshift of interest. We 
construct Monte-Carlo realizations of "merger trees" using 



the method of somerville & Kolatt (1998), which was tested 



agains t merger trees extracte d from dissipationless simula- 
tions (Somerville et al. 2000a). Each branch in the tree rep- 



resents a halo merging event, and the trees are truncated 
when the circular velocity of the progenitor halo becomes 
smaller than 40 km s _1 . We assume that gas in halos smaller 
than this effective mass resolution is photoionized and can- 
not cool or form stars (see SP). We construct merger histo- 
ries for a grid of halos with circular velocities ranging from 



40 km s _1 to 1500 km s _1 , and weight the results using the 
appropriate Press-Schechter probability for the appropriate 
halo mass and redshift. 

When a halo collapses or merges with a larger halo, we 
assume that the associated gas is shock heated to the virial 
temperature of the halo. This gas then radiates energy and 
cools. The cooling rate depends on the density, metallicity, 
and temperature of the gas. Cold gas is turned into stars us- 
ing a simple recipe, and supernovae energy reheats the cold 
gas according to another recipe. To model the star formation 
rate, we use an expression of the general form 

m, = , (1) 

where m* is the star formation rate, m co id is the total mass 
of cold gas in the disc, and cold gas is converted into stars 
with a time scale t*. In principle this time scale could be 
a function of the circular velocity or dynamical time of the 
disc or other variables. In the simplest version of this recipe, 
t» = t° = constant; i.e. the efficiency of conversion of cold 
gas into stars m*/m co id is independent of the other proper- 
ties of the galaxy. This recipe was one of those investigated 
in SP (called SFR-C in the terminology of SP; here called 
"constant efficiency quiescent"). We also consider a recipe 
in which r* oc id yn , where id yn is the dynamical time of 
the galaxy. Because the density of the collapsed haloes is 
higher at high redshift, the typical dynamical times become 
smaller and the conversion of gas into stars becomes faster. 
This recipe is similar to the one usually used by the Munich 
group (called SFR-M in SP), and here we call it "accelerated 
quiescent." 

Supernovae feedback is modelled using the disc-halo 
model introduced in SP. The rate of reheating of cold gas is 
given by 



ESN 



(2) 



where esn is the rate at which energy is injected into the cold 
gas by supernovae, and Vg SC is the average escape velocity of 
the disc or halo (the reheating rate and ejected gas mass is 
calculated seperately for each of these components). 

Each new generation of stars produces a fixed yield of 
metals, which are immediately deposited back into in the 
cold gas. The metals may then be mixed with the hot gas in 
the halo, or ejected from the halo by supernovae feedback. 
This in turn affects the cooling, which becomes more efficient 
as the metallicity of the gas increases. Unlike our work in SP, 
where we used a fixed metallicity for the hot gas in calcu- 
lating the cooling function, here we use the modelled value 
of the hot gas metallicity and the me tallicity- dependent ra- 



diativ e cooling curves tabulated by Sutherland & Dopita 



(199$ . 



Another new feature (not present in the models of SP), 
is the treatment of gas and metals that are ejected from 
the halo. In SP this material was never re-incorporated in 
any halo. Here, we assume that the material is distributed 
outside of the halo with a continuation of the isothermal 
r~ 2 profile that we assume inside the halo, and such that if 
the total mass of the halo were to double, all of the mate- 
rial would be re-incorporated in the halo. This material falls 
in gradually (according to a spherical infall model) as the 
virial radius of the halo increases due to the falling back- 
ground density of the Universe. This is similar to the recipe 
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used by CAFNZ and Kauffmann et al. (1998), in which all 



of the hot gas ejected from the halo is re-incorporated (all 
at once) when the halo doubles in mass, but it removes the 
discontinuous behaviour caused by the sudden infall of a 
large amount of gas. The assumption that the ejected gas 
should return when the halo mass doubles is just an arbi- 
trary relic of the block model used by CAFNZ, in which the 
halo masses always grew by factors of two in each branch 
of the merging tree; the actual timescale for return of the 
ejected gas is quite uncertain. 

We now also account for the mass that is returned to the 
ISM by young stars. For each quantity of mass that is turned 
into stars, a fraction 7? is returned to the cold gas reservoir. 
Here we assume R = 0.14, which results from assuming that 
all mass from stars with masses greater than 8 Mq is re- 
cycled, using a Salpeter IMF with an upper mass cutoff of 
100 Mq. Had we also included mass loss from smaller-mass 
stars, the value of R would be larger, but these estimates 
are more uncertain so we neglect this contribution. 

2.2 Mprgprs anH Starhnrsts 



lark matter haloes merge, the galaxies contained 
hem survive and may merge on a longer time scale, 
halo merger event, the central galaxy of the largest 
cni|tor halo becomes the new central galaxy, and all 



When 
within 
After a 
prog. 

other g ' alaxica become "3atcllitc3." — Satellite galaxies may 

merge with the central galaxy on a dynamical friction time 
scale, or with each other on approximately a mean free path 
time scale. However, if the relative velocities of the satellites 
are large compared to their internal velocities, they will not 
experience a binding merger. Thus the satellite merger rate 
decreases in clusters. The expressions for the satellite-central 
and satellite-satellite merger rates are given in SP. 

We assume that every galaxy-galaxy merger triggers a 
starburst. Our treatment is based on a simple parameteriza- 
tion of the results of N-body simulations with gas dyn amics 
whic h we now summarize briefly. Mihos & Hernquist (1994 



199f) have simulated galaxy-galaxy mergers using a high 



resolution N-body/smoothed particle hydrodynamics code 
(TREESPH) with star formation modelled according to a 
Schmidt law (psfr oc p™ as with n = 1.5). In Mihos & Hern- 



quist (1996), mergers between equal-mass galaxies (major 
mergers) were simulated, and it was found that 65-85 per- 
cent of the total gas supply (in both galaxies) was converted 
into stars over a time scale of 50-150 Myr. These results were 
not very sensitive to morphology or the orbital geometry. 
The case of highly uneq ual-mass mergers was explored 



Mihos & Hernquist (1994). The case simulated represents 



a Milky Way-sized disc galaxy accreting a satellite that is 
one tenth of its mass. The non-axisymmetric mode gener- 
ated by the accretion of the satellite causes a large fraction 
of the gas to collapse into the central region of the galaxy, 
fueling a strong starburst. In this case about 50 percent of 
the original gas supply is consumed in the starburst, which 
lasts for about 60 Myr. However, the results are much more 
sensitive to the morphological structure of the galaxies than 
the major-merger case. If the larger galaxy has a bulge (the 
case simulated has a bulge to disc mass ratio of 1:3), the 
bulge seems to stabilize the disc against strong radial gas 
flows, leading to a much weaker starburst event (o nly about 
5 percent of the total gas supply is consumed). Mihos & 
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Figure 1. 

of the bar 
are scaled 



and Mihos 



The efficiency of star formation in bursts as a function 
'onic mass ratio of the mer ging galaxies. The results 
f o match the calcul ations of Mihos & Hernquist (1994) 
fc Hernquist (199^ ). When a bulge is present, the burst 



efficiency c rops off more quickly at lower mass ratios. 



Hernquist (1994) note that this implies that the importance 
of bursts will decrease at low redshift as galaxies develop 
bulges, even if the merger rate remains constant. 

These results inform our treatment of collisional star- 
bursts. When any two galaxies merge, the "burst" mode of 
star formation is turned on. The star formation rate due to 
the burst is modelled as a Gaussian function of time with a 
width of (Tburst- The burst is switched off after a time 4o"burst 
has elapsed. The burst model has two parameters, the time 
scale and the efficiency of the burst. The efficiency eburst 
is defined as the fraction of the cold gas reservoir (of both 
galaxies combined) that is turned into stars over the entire 
duration of the burst. We assume that the efficiency is a 
power-law function of the mass-ratio of the merging galaxy 
pair: 



Gburst 



/con 



[msmaii/rribig 



with the two paramet ers /, 



consume 



(3) 

and Qburst chosen to re- 



merger, 

(/coj 



and 



produce the results of Mihos fc Hernquist (1996 ) for a 1:1 



Mihos & Hernquist (1994) for a 1:10 merger 



onsumc — 0.75 and ftburst = 0.18). When a bulge of at least 
a third of the disc mass is present, bursts are suppre ssed in 
minor mergers (based on Mihos fc Hernquist 199^ ). Fig.l 
shows the resulting scaling of the burst efficiency with the 
mass ratio of the merging pair for the bulge and no-bulge 
cases. We assume that the burst timescale e>"burst is propor- 
tional to the dynamical time of the larger of the two discs. 
The quiescent mode of star formation continues uninter- 
rupted according to Eqn. ^ above, although its contribution 
is generally insignificant in comparison to the burst. 

Each merger is classified as "major" or "minor" ac- 
cording to whether the ratio of the smaller to the larger 
of the galaxies' baryonic masses is greater than or less than 
the value of the parameter /bulge ~ 0.25. Major mergers 
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have mass ratios greater than /bui gc , and the bulge and disc 
stars of both galaxies plus all new stars formed in the burst 
are placed in a bulge component. Minor mergers have mass 
ratios less than /bulge, and the pre-existing stars from the 
smaller galaxy are placed in the disc component of the post- 
merger galaxy while all newly formed stars again join the 
bulge. This is motivated by the N-body simulations of simi- 
lar satellite accretion events by Walker, Mihos, & Hernquist 
(1996^71 which show that 90 percent of the satellite mass is 
stripped and distributed in the disc of the larger galaxy and 
only the densest part of the satellite core ends up at the 
centre. 

A possible problem in scaling from the z = simula- 
tions is that galactic discs may be quite different at high 
redshift. Galaxies then are smaller and denser and proba- 
bly have higher ratios of gas to stars. This may significantly 
affect their susceptibility to star-forming instabilities. We 
are currently investigating this using a more extensive set 
of N-body hydrodynamic simulations similar to those of Mi- 
hos & Hernquist, but with the initial galaxy prop erties cho- 
sen to be rep resentative of high-redshift galaxies ( somervilk 
et al. 2(l00b). Similarly, the modelling of mergers, infall and 
tidal stripping of sub-haloes should be refined and tested 
against high-resolution dissipationless cosmolo gical Simula- 



rcdshiT 0.4 



2 3 4 5 10 



tion s. W ork on these issues is also in progress (Kolatt et al 
200(f~| p 

Fig. H shows the total star formation rate for the 
largest progenitor of the central galaxy in a Milky Way- 
sized halo (V c = 220 km s" 1 ) and in a group-sized halo 
(V c — 500 km s _1 ), both at z — today. The star formation 
rate is shown in models with: (1) no starbursts, (2) bursts in 
major mergers only, and (3) bursts in both major and minor 
mergers. All three models contain quiescent star formation 
using the constant efficiency recipe. From this figure we can 
see that minor mergers are important even if their efficiency 
is low because they are much more common than major 
mergers. Comparing the group-sized halo with the Milky 
Way halo shows how mergers shut down in high-velocity- 
dispersion environments. This is because we have assumed 
that the relative velocity of colliding objects must be small 
compared to their internal velocities in order for objects to 
merge. The burst models discussed in the remainder of this 
paper correspond to model (3) unless specified otherwise. 



2.3 Stellar Population Synthesis 

With the full star formation history of a galaxy in place, we 
estimate its luminosity in any desired passband using the 
stellar population synthesis models of Bruzual & Chariot 
(GISSEL). In this paper, our standard treatment uses the 
1998 version of the solar-metallicity models with a Salpeter 
IMF. Modelling the young, massive, sub-solar metallicity 
stars that might be found in LBGs is probabl y quite un- 
certain (see [Chariot, Worthey, fc Bressan 1996 ). However, 
it is encouraging that several groups have produced models 
that agree quite well in the UV and optica l part of the SEP 



(see Devriendt, Guiderdoni, & Sadat 1999 



and Section 



5.2) 



The mass-to-light ratio in the UV is relatively insensitive to 
metallicity for young stars, but it is quite sensitive to the 
IMF (primarily to the ratio of high-mass to low-mass stars) . 
We include only starlight and neglect the contribution from 
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J -2 
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redshift 0.4 1 2 3 45 10 



quiescent 



pro:o-"group' 



major bursts 
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5 10 

lookback time (Gyr) 

Figure 2. The total star formation rate for the largest progenitor 
of the central galaxy within a halo with a present-day circular ve- 
locity of 220 or 500 km s — 1 . All three models contain constant ef- 
ficiency quiescent star formation. The top panel shows the model 
with no bursts, the middle panel shows the model with bursts in 
major mergers only, and the bottom panel shows the models with 
bursts in major and minor mergers. 
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nebular emission. See Section for a detailed discussion of 
uncertainties in the stellar population synthesis. 



2.4 Dust 

Because most of the observations of high-redshift galaxies 
are obtained in the rest-UV, the effects of dust are likely 
to be important. Just how important has been a matter 
of debate ever since the discovery of the LBG population. 
Pettini et al. (1997a) estimated an extinction at ~ 1500 A 
of a factor of ~ 3, based on the ratio of emission lines to 



the continuum in a f ew objects. Meurer et al. (1997 ) and 
Sawicki & Yee (1998) estimated much larger factors of 15- 



20. Recently, there seems to have be en convergence to an 
intermediate value of a factor of ~ 5 (Meurer, Heckman, & 



Calzetti 1999 



Steidel et al. 199S| ). This work makes use of a 



correlation between the FIR excess (a reliable observational 
measure of bolometric extinction) and the far-UV spectral 
slope in nearby starburst galaxies. If the same correlation is 
then assumed to hold in high-redshift galaxies, the measured 
UV spectral slopes of LBGs indicate that there is an average 
extinction of a factor of 4.7 in t he relatively bright (TZ < 
25.5) LBG populatio n studied by fcteidcl et al. ( 1999b - The 
same measurements (Meurer, Heckman, & Calzetti 199£ ) 
indicate that, in LBGs as in local starburst galaxies, the 
most UV-luminous (and rapid ly star forming) objects are 
the most heavily extinguished (Wang & Heckman 1996). 

To estimate the effects of dust on our model galaxies, we 
use a very simple parame terization of the empirical results of 
Wang & Heckman (1996) for nearby starburst galaxies. The 



optical depth for a galaxy with an intrinsic (unextinguished) 
UV luminosity Luv,i is given by: 



tuv = Tuv,t 



Li 



~*UV,i 



(4) 



This is identical to the recipe used in SP except that we now 
normalize the recipe in the UV rather than the B band. We 
take Luv.* to equal the observ ed value of L* in the z — 3 
sample of Steidel et al. (1999] ), i.e., the luminosity corre- 
sponding to niAB = 24.48 in the appropriate cosmol ogy. We 
use the sam e value of L, at redshift z = 4, which as Steidel 
et al. (1999[ ) note is consistent with the z = 4 data. We take 
P = 0.5 as in Wang & Heckman (1996). We then assign to 
each galaxy a random inclination and compute the actual 
extinction using a standard slab model (see SP). We adjust 
the parameter controlling the face-on optical depth, tuv,*, 
to obtain an average extinction correction at 1500 A of a 
factor of ~ 5 for z ~ 3 galaxies with luminosities typical 
of the Steidel et al. sample, consistent with the estimates 
described above. The value tuv,* = 1-75 gives good results. 
To extend the extinction correction to other wavebands, we 



assume a Calzetti attenuation curve (Calzetti 1997) 



2.5 Model Parameters 

Throughout this paper, unless stated otherwise, we use the 
fashionable ACDM cosmology with Qo = 0.3, Qa = 0.7, 
h = # /(100km/s/Mpc) = 0.7, a 8 = 1.0. These values are 
consiste nt with a great deal of current data (for a recent re- 
view see Primack 200C ) . As in SP, the main free parameters 



related to galaxy formation — those describing the quies- 
cent star formation efficiency (r«), the supernovae feedback 
efficiency (egjj), and the metallicity yield (y) — are set by re- 
quiring an average present-day "reference galaxy" (the cen- 
tral galaxy in a halo with a circular velocity of 220 km s _1 ) 
to have (a) an I-band magnitude Mi — 5 log h = —21.8, (b) a 
cold gas mass m co id — 10 10 h~ 2 M©, and (c) a stellar metal- 
licity of about solar. Requirement (a) adjusts the zero-point 
of the I-band Tully-Fisher relation to agree with observa- 
tions, while requirement (b) fixes the cold gas content of the 
"reference galaxy" to agree with observed Hi masses, allow- 
ing for a contribution from molecular hydrogen (see SP). 
The values of the free parameters used here are similar to 
those used in SP. 

In some previous work (e.g., BCFL), it has been as- 
sumed that only a fraction f* nm < 1 of the stars formed in 
the models are luminous, the remainder being in the form 
of brown dwarfs or other non-luminous baryonic material. 
We find that in order to get our reference galaxy to be 
bright enough today with our assumed Salpeter IMF, wc 
need ff w ~ 1. 

Note that the parameters have been set entirely by com- 
paring to a subset of present-epoch observations. We have 
shown in SP that the models then reproduce many impor- 
tant observed features of nearby galaxies. Now, leaving all 
free parameters fixed, we can consider the predictions of the 
same models for the high-redshift Universe. 



3 MODEL RESULTS 

Because extensive spectroscopic follow-up work has now ver- 
ified the effectiveness of the Lyman-break or "drop-out" 
technique for finding z > 2 galaxies, we do not attempt 
to mimic the colour selection process used to identify real 
Lyman-break objects but rather assume that all galaxies 
brighter than a certain limiting magnitude in our models 
would be in fact be selected as LBGs. The modelling work 
of BCFL has shown that using the same colour selection 
criteria as the observations picks out model galaxies in the 
expected redshift range. We do not include the effect of ab- 
sorption by intervening cold gas clouds because, although 
this effect is very important shortwards of the Lyman limit, 
for redshifts less than z ~ 4 it does not much affect the 
spectral energy distribution (SED) in the wavelength range 
relevant to our study. We have calculated magnitudes using 
the filter response functions of the WFPC2 filters F606W 
an d F814W (Veoe and Isi-i), as well as the 1Z and I filters 
of Steidel & Hamilton (1992) used in the ground based ob- 



servations and kindly provided to us in electronic form by 
K. Adelberger. All m agnitudes are given in the AB system 
( |Oke fc Gunn 1983; ). The U 60 6 filter and the K filter cor- 
respond to a mean rest wavelength of 1600 A at redshifts 
2.75 and 3.04 respectively. The 7si4 and Steidel & Hamilton 
I filters correspond to the same rest wavelength at redshifts 
of about 4.00 and 4.13 respectively. 

We now consider three different recipes for star forma- 
tion. The model referred to as the "collisional starburst" 
model includes quiescent star formation using the "constant 
efficiency" (CE) recipe, plus star formation in bursts as de- 
scribed in the previous section. The "constant efficiency qui- 
escent" model has no burst mode, and all star formation 
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LBGs with spectrosc opic redshifts, we ha ve used the values 
given in Table 3 of |Steidel et al. (1999|). The la field-to- 
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Figure 3. The comoving number density of galaxies brighter than 
25.5 at rest f500 A. Thin lines show the three models with no 
correction for dust extinction, and bold lines show the models 
with du st extinction included using the recipe described in Sec- 
tion 2.4 . The filled squares indicate the comoving number density 
of LB Gs with spectroscop ic redshifts from the ground-based sam- 
ple of Btcidcl et al. (1999| ). Four- pointed stars indicat e the values 
derived from LBGs in the HDF ( Pozzctti ct al. 199 j ). The stars 
and triangles are from the HDF-N and HDF-S gal axies with pho- 
tom etric redshifts from the Stonybrook catalogs ( [Lanzetta et al 
1999"" 



l etric rcdshilts trom the btonyprook catalogs ( |Lanzetta et al. 
offJl icn ct al. 199S| ). The observations have not been corrected 



for extinction. Note that the constant efficiency quiescent model 
predicts a strong decrease in the number density of bright galax- 
ies with redshift, in contrast to the other two models, which show 
a very gentle decline from z ~ 2 to z ~ 6. 



is modelled using the constant efficiency recipe. The third 
model is also quiescent, but star formation efficiency varies 
inversely as the dynamical time of the galaxy. We refer to 
this third model as "accelerated quiescent" because the star 
formation rate for a given cold gas mass becomes higher at 
high redshift due to the increasing density of the haloes. 
Note that the supernovae feedback recipe has been kept the 
same and all other free parameters have been left fixed to 
the same values for the three models. In SP, we showed 
that, by adjusting only the normalization of the star for- 
mation efficiency, r» , it is possible to make all three models 
to match observed properties of local galaxies such as the B 
and K band luminosity functions, the Tully-Fisher relation, 
colours, gas and metal contents, etc. 



3.1 Comoving Number Density 

We first investigate the model predictions for the number 
density of bright objects as a function of redshift. Fig. ^ 
shows the comoving number density of galaxies brighter 
than an apparent AB magnitude limit of mAs(1500) = 25.5 
at rest 1500 A over the redshift range 2 < z < 6. The co- 
moving number densities for the observations have been cal- 
culated using the appropriate geometry for the ACDM cos- 
mology used in our models. For the ground-based sample of 



field variance of this sample is about 12 percent; if plotted 
on Fig. ^, the error bar would be smaller than the sym- 
bol. The comoving number density of LBGs in the HDF has 
been calcu lated using the Veoe counts of U-dropouts from 
Table 1 of Pozzetti et al. (1998), and the 7~8i4 counts of B- 



dropouts helpfully provided to us by L. Pozzetti. Error-bars 
on these points represent Poisson errors only. We also show 
the number densities computed from the Stonybrook cat- 
alog of photometric r edshifts from th e HDF-N and HDF-S 
( Lanzetta et al. 1999| ; Chen et al. 1999 ) , generously provided 
to us by H.-W. Chen and K. Lanzetta. Once again the error 
bars are simply la Poisson and do not account for photo- 
metric or redshift errors. Note that, unlike the point derived 
from galaxi es identified in the H DF using the Lyman-break 
technique (Pozzetti et al. 199c), the values from the pho- 
tometric redshift catalogs at z ~ 3 agree with the ground- 
based estimate within the errors. This suggests that the full 
photometric redshift technique, including the observed near- 
IR bands, finds high-redshift galaxies that are missed by the 
Lyman-break technique. The difference between the HDF-N 
and HDF-S (stars and triangles) gives an indication of the 
field-to-field variance characteristic of the HDF volume. 

Fig. |3] demonstrates an important and robust predic- 
tion of the models — the constant efficiency quiescent mod- 
els predict a strong decline in the number density of bright 
galaxies with increasing redshift. This is a generic feature of 
any hierarchial model in which the efficiency of the conver- 
sion of cold gas into stars remains constant with redshift. 
This is because the number density of the massive haloes 
that host bright galax ies in the quiescent model drops off 
sharply with redshift. Kolatt et al. (199£) found a similar 
result based on iV-body simulations. The net result is that 
the CE quiescent model without dust extinction produces an 
acceptable match to the number of bright galaxies at z ~ 3 
but badly underpredicts the numbers at redshifts z > 4. 
The shortfall is unacceptable at all redshifts when dust is 
included. 

In contrast, the burst model contains collisional star- 
bursts that cause smaller-mass objects to become bright 
enough to exceed the detection limit for a brief amount of 
time. Because the collision rate is relatively constant over 
the redshift range 2 < z < 6 and there is an ample supply 
of cold gas, the comoving number density of visible galax- 
ies remains almost constant over this range. The collisional 
starburst model actually overpredicts the number of bright 
galaxies, but adding dust provides excellent agreement at all 
redshifts. Finally, the accelerated quiescent model produces 
nearly identical results to the collisional starburst model, 
but for a very different reason. Instead of being temporarily 
brightened by bursts, each galaxy has a constant low mass- 
to- light ratio because of the accelerated star formation rate. 
This model also does a good job of matching the number of 
moderately bright galaxies at all redshifts. 

The observations at 2 <C 2 <J 4.5, especially those base d 
on the large ground-based sample of steidel et al. (1999), 



provide a secure lower limit on the number density of bright 
galaxies. The situation at higher redshift is more uncertain, 
but the differences between the constant efficiency quiescent 
model and the other two are quite dramatic. According to 
the CE quiescent model, the probability of finding even one 
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Figure 4. The luminosity function at z = 3 in the observed 
IZab band (rest ~ 1600 A). Symbols show the composite lumi- 
nosity function fro m ground-based observations (filled squares 
Steidel et al. 1999 and the HDF (filled circles; |Pozzetti et al 



1998 |). (j )pen squares show the luminosity function derived from 
the ground-based observations after correcting for the effe cts of 



dust ex tinction using the method described in the text (Adcl- 
berger <S ; Steidel 2000,). Lines indicate the results of the collisional 



starburst model without (solid) and with (dashed) corrections for 
dust extinction. 



bright galaxy at z J> 5 in a volume comparable to the HDF is 
only one in 10 6 , whereas both the collisional starburst model 
and the accelerated quiescent model predict on the order of 
a few galaxies per HDF at 4.5 < z < 6, even with dust. A 



been discovered recently (Dey et al. 1996; Weymann et al. 


1998]; |Spinrad et al. 1998 


; Chen, Lanzetta, & Pascarelle 19991 


van Breugel et al. 1999; 


Hu, McMahon, & Cowic 1999), but, 



as most of these objects were found serendipitously, it is 
difficult to estimate their number density. Such candidates 
are also being detected in fair numbers in photometric red- 
shift catalogs, but these detections are also uncertain, and 
further spectros copic work will be requ ired in order to con- 
firm them (see Stern & Spinrad 1999 for a recent review 
on search techniques for very high-redshift galaxies). If the 
current preliminary detections of very high-redshift galax- 
ies are confirmed, they will essentially rule out the constant 
efficiency quiescent model. 



3.2 Luminosity Function 

The luminosity function and its evolution with redshift pro- 
vide additional tests of models. The bright end (jtiab < 
25.5) o f the luminosity function at ~ 1500 A for the ob- 
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Figure 5. Same as Fig. M, for the quiescent models. 



resultant composite luminosity function at z — 3 is shown in 
Figs. ^ and |E| The ground-based data have been corrected 
for incompleteness (see Steidel et al. 1999). As usual, the 



observations have been converted to our ACDM cosmology. 

Fig. |^ shows that the luminosity function of the colli- 
sional starburst model with dust has a steep, Schechter-type 
drop-off at bright magnitudes, whereas the intrinsic func- 
tion without dust is much flatter. Thus in this picture, the 
knee seen in the observations is purely an artifact of dif- 
ferential dust extinction. It is gratifying that our intrinsic 
luminosity function is very similar in shape to t he observed 
dust-corrected function derived from the data (Adelberger 



served ^BGs at z ~ 3 can be determined fairly accurately fc Steidel 200C ), for which the corrected luminosity of each 



from the relatively large ground-based spectroscopic sam- 
ple. Th e HDF can then be used to give some indication of 



galaxy has been derived using the measure d far-UV slope 



and the dust-slope cor relation discussed in Meurer, Heck 



the fair t end slope, although there is pro bably s ome incom- 



pleteness faintwards of about m = 26 ( Dickinson 1996 ) . The 



man, fc Calzetti (1999) 



Fig. H shows similar plots for the two quiescent models. 
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Figure 6. Same as Fig. for the I band at z = 4 (rest ~ 1600 
A). Observed data arc shown as solid symbols (see key on Fig. ^|). 
The dotted line shows the Sche chter fit to the 2 = 4 luminosity 
function of Steidel et al. (1999), which is somewhat higher than 



the HDF results. The data may be affected by incompleteness 
below ~26 magnitude. 



The constant efficiency model without dust (bottom panel) 
has a steeply decreasing luminosity function at the bright 
end, reflecting the steepness of the halo mass function at 
these large masses. When dust extinction is added, the CE 
quiescent model becomes even steeper and dramatically un- 
derpredicts the number of bright galaxies, consistent with 
its underprediction of count data in Fig. ^. The situation 
with the accelerated quiescent model is rather different. Al- 
though that model predicts about the same total number of 
bright (m^B (1500) < 25.5) galaxies as the collisional star- 
burst model in Fig. ^, the shape of its luminosity function is 
different at the very brightest magnitudes. The burst model 
shows a tail at bright magnitudes, whereas the accelerated 
quiescent model even without dust drops off steeply, like a 
Schechter function. With dust included, the accelerated qui- 
escent model falls even more steeply and fails to match the 
brightest galaxies. This, again, is due to the tight link be- 
tween halo masses and galaxy luminosities in all quiescent 
models. 

Figs. |^ and (?] show analogous comparisons for luminos- 
ity functions at z ~ 4. The collisional starburst model fits 
the luminosity function for bright galaxies whereas both qui- 
escent models fall short, especially when dust is added. All 
models predict too many faint galaxies. This may be due to 
the fact that, with the parameters we have chosen, our dust 
model predicts very little extinction in galaxies fainter than 
m ~ 26. In correcting the models for dust, we have simply 
assumed the same parameters as at z ~ 3. Note also that 
the observed luminosity function faintwards of m ~ 25 (the 
magnitude reached by the ground-based survey) is quite un- 
certain and possibly incomplete. Both effects would tend to 
make the models deviate in this way. 

In summary, the CE quiescent models predict a strong 




Figure 7. Same as Fig. for the quiescent models. 

evolution of L* with redshift, wheras the burst model and 
the accelerated quiescent model predict very little evolution 
in L* over the redshift range 3 < z < 5. The observed and 
dust- corrected luminosity functions of the collisional star- 
burst model are in good agreement with all data at both 
redshifts. However, systematic errors in the dust correction 
could be large. A more direct test of the number densities 
of very bright galaxies (which correspond to objects with 
extremely high star formation rates) will come from future 
sub-millimeter observations. 

3.3 The Star Formation History of the Universe 

The star formation history of the Universe, i.e., the global 
rate of star formation as a function of redshift, is a crucial 
test of theories of galaxy formation and cosmology. There 
are many observational tracers of the star formation his- 
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Figure 8. The star formation history of the Universe ("Madau 
diagram") for our ACDM cosmology. Symbols indicate observed 
estimates of the star formation rate density at various redshifts; 
sources of data are listed in Table A2. All data points have been 
correcte d for incompleteness but not for dust extinction. Lines 



show pr ;dictions of the three models. 



tory, each with their own selection effects. In this section we 
investigate the star formation history traced by the lumi- 
nosity density, age constraints on stars in nearby galaxies, 
cold gas detected in quasar absorption systems, and metals. 

3.3.1 Luminosity Density 

The star formation history of th e Universe, as repr esented by 
the popular "Madau diagram" (Madau et al. 1996), is an im- 
portant , signatm3Mjf_ai^ 



ever, ai accurate determination of this diagram is far from 
straightforward: the observations must be converted from 
luminosity to star formation rate and corrected for incom- 
pleteness and dust extinction (see Appendix). As new ob- 
servations have been added at various redshifts, the diagram 
has undergone a continuing metamorphosis from its original 
form. Fig. H shows the star formation history obtained from 
our burst and quiescent models, along with a compilation of 
recent observations, where the observations have been cor- 
rected for incompleteness but not for dust extinction. The 
observations have been converted to our ACDM cosmology 
as described in the Appendix. Note that the constant ef- 
ficiency quiescent model shows the characteristic peak at 
z ~ 1.5 and the relatively steep decline at higher redshifts 
t hat were notable f eatures of the original Madau diagram 
(Madau et al. 1996). However, the data points in that dia- 



gram at z = 2.75 and z — 4 were based on the very small 
HDF s urvey and relied on the Lyman-break technique to 



select galaxies in the desired redshift range. As Madau et 
al. correctly emphasized, these points should be regarded as 
lower limits. More recent observations based on much larger, 



Figure 9. The new Madau diagram, corrected for dust extinction. 
Symbols indicate observational estimates of the star formation 
rate density at various redshifts. Triang les are derived from Ha 
(Gronwall 199£ ; Trcsse & M addox 1998). and s q uares arc derived 
from rest-UV l u minosities (Trever et al. 1998 |Cowie, Songaila.' 



& Barger 199E; Steidel et al. 1999). These data have been cor 



rected for incompleteness and dust extinction (see Appendix for 
details). The upper-limit symbols show the data with a "max- 
imal dust correction," in which a fixed factor of five correction 
has been applied to all galaxies. If dust extinction is luminosity 
dependent, this probably overestimates the true extinction cor- 
rection (see text). Diamonds indicat e results derived fr om far-IR 
sources at z = 0.7 observed by ISO (Flores et al. 1999) and sub- 
mm sour ces, believed to be at high redshift (~ 3), observed by 



SCUBA (Hughes et al. 1998). Lines show the three usual models. 



ground-bas ed surveys with spectroscopic redshifts (Steidel 



et al. 1999 ) show a much shallower evolution between z 



and 2 ~ 4. 

Another important recent observational development is 
a more convincing estimate of the effects of dust extinction 
on the values derived from luminosities in the far-UV. This 
has also produced a si gnificant modificatio n of the diagram, 
as recently noted by Steidel et al. (1999). Fig. ^| shows a 
subset of the observations with a correction for dust extinc- 
tion based on the observational results (see Appendix for 
details). The CE quiescent model is now seen to be highly 
inconsistent with the revised observational estimates at high 
redshift. The collisional starburst model shows a gentler rise 
and a broad peak around z ~ 4, with near-constant star for- 
mation density out to a redshift of ~ 6. This appears very 
consistent with the recent optical estimates after correction 
for dust extinction (see Appendix); it is also consistent with 
the luminosity densities implied by the far-IR and s ub- mm 
sources detected at intermedia te and high redshift (Flores 



et al. 1999; Hughes et al. 1998). The redshift dependence of 
the star formation rate density in the accelerated quiescent 
model is very similar to the collisional starburst model, al- 
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Figure 10. The star formation rate in the collisional starburst 
model is shown broken down into the contribution from quiescent 
star formation and bursts. The bursting mode dominates at high 
redshifts, and the quiescent mode dominates at low redshifts. The 
cross-over occurs at a redshift of about 0.8. 
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Figure 11. Star formation in the bursting mode is shown bro- 
ken down into contributions from different kinds of merger events. 
The dotted and dashed lines show contributions from major and 
minor mergers (respectively) between a satellite galaxy and a cen- 
tral galaxy. The long-dashed and dot-dashed lines show the con- 
tribution from major and minor mergers between two satellite 
galaxies. The dominant contribution at all redshifts comes from 
minor mergers onto central galaxies. 



Figure 12. The density in stars, in units of the critical den- 
sity (f2*). The observational estimate at z = is shown by the 
symbol with the errorbar. The observational constraint at high 
redshift from "fossil evidence" (see text) is shown by the shaded 
box. The collisional starburst model and the accelerated quies- 
cent model produce a significant fraction of their stars at high 
redshift, wheras the constant efficiency quiescent model produces 
stars too late. 

though the overall star formation rate density is a bit lower 
than the data. 

Fig. ^ shows the contribution to the star formation rate 
density from the burst and quiescent modes separately for 
our collisional starburst model. The bursting mode domi- 
nates at high redshift but declines more steeply than the 
quiescent mode as redshift decreases. A crossover therefore 
occurs, at z ~ 0.8, and quiescent star formation dominates 
the low-redshift Universe. 

Fig. |ll] shows the contribution to star formation in the 
burst mode from different kinds of merger events. We show 
the contribution from major and minor mergers between 
satellite and central galaxies, and major and minor mergers 
between two satellite galaxies. Note that in previous semi- 
analytic models, starbursts were included only in the first 
case, i.e., satellite-central/major mergers. We can see from 
the figure that the contribution from other kinds of merg- 
ers, particularly satellite-central/minor, is quite important. 
The relatively small contribution of satellite-satellite merg- 
ers may be an artifact of our rather questionable simplifying 
assumption that all new cold gas is accreted onto the central 
galaxy; this assumption tends to deprive small satellites of 
the fuel needed to produce bright starbursts, but may not 
be realistic. 

3.3.2 Fossil Evidence 

Another way of constraining the star formation history of 
the Universe is from the "fossil evidence" contained in the 
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ages of stars in present-day galaxies. Renzini (1996) has 



argued that constraints on the ages of early-type galaxies 
in clusters from the small observed scatter of Fundamental 
Plane and colour-magnitude relations, combined with the 
fraction of today's stellar mass contained in early-type sys- 
tems (bulges and ellipticals), can be used to deduce that one 
third of the stars we see t oday must have formed at z > 3. 
Van Dokkum et al. (1998) find that the formation redshift 



of early-type galaxies in clusters may be relaxed to z J> 2 
in cosmologies with a large cosmological constant (such as 
our ACDM cosmology) because of the longer time that has 
ellapsed in such models. Fig. [l^ shows the density of stars 
in units of the critical density for the burst and quiescent 
models. The point with the lar ge error bar at redshift zero 



is from the "baryon budget" of Fukugita, Hogan, fc Peebles 
( 1998^7] All three models agree with this estimate within the 
errors. The shaded patch shows the constraints on f2„ at 
high redshift from a Renzini-like argument, where we have 
simply divided the stellar baryon budget at z = (and its 
error bar) by a factor of three. The CE quiescent model 
does not produce enough stars at high redshift, but the col- 
lisional starburst and accelerated quiescent models meet this 
constraint. The CE quiescent model produces one-third of 
the stars by z = 1.5, and half of the stars by z = 0.9. The 
collisional starburst and accelerated quiescent models pro- 
duce one-third of their stars by z = 2.2, and half of their 
stars by z — 1.5. 



3.3.3 Cold Gas Density 

Observations of quasar damped Lyman-a systems (DLAS) 
provide an estimate of the Hi and metal content of the Uni- 
verse from z ~ 0.7 to z ~ 4. Fig. |l^ compares the models 
to an estimate of the fraction of the crit ical density in the 



form of cold gas from the observations of storrie-Lombardi 



McMahon, & Irwin (1996) (the data converted to the ACDM 



cosmology was kindly provided by R. McMahon and C. Per- 
oux). The observations shown should be considered a lower 
limit on the total mass of cold gas for several reasons. Some 
dusty DLAS may also be missed because their background 
quasars would be too dimmed to be included in an o ptically 
selected, magnitude limited sample (Pei & Fall 1995). In ad- 
dition a significant amount of mass could be in the form of 
lower column density Hi clouds or in the form of molecular 
or ionized hydrogen. 

We can see from Fig. [l3] that the CE quiescent model 
(dashed line) overproduces the amount of cold gas by a fac- 
tor that is perhaps a bit large to be explained by these ef- 
fects, whereas the collisional starburst model (solid line) is 
more consistent with the data. The accelerated quiescent 
model (dotted line) is roughly la lower than the data at 
z > 2, and this may be a concern since all of the effects 
mentioned above will tend to increase this discrepancy. Note 
that the constraint provided by the DLAS is complemen- 
tary to that provided by observations of LBGs. We could 
decrease the consumption of cold gas by decreasing the effi- 
ciency of star formation r« , but this would in turn produce 
fewer bright LBGs. 
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Figure 13. The evolution of f2 co id as a function of redshift. Data 
points show the matter density in Ht from ob servations of DLAS 
(Storric-Lombardi, McMahon, & Ir win 1996). The d ata point at 



2 = is from local Hj observations ( Zwaan et al. 1997 ). The quan- 
tity plotted for the models (line types are as in Fig. |12|) is generic 
"cold gas", some of which may be in the form of molecular or 
ionized hydrogen, or may be in dust extinguished or low-column 
density systems that would not be detected as damped systems. 
Thus the model lines are upper limits on the quantity that is 
actually determined from the observations (f2m in DLAS) and 
should lie above the data points. Line types are the same as in 
figure p| 



3.34 Metals 

The metal content of the Universe is yet another tracer of the 
star formation history. The mean metallicity of the DLAS 
can be estimated from Znn and Cm absorption lines, which 
are not much depleted by dust (Pettini et al. 19971). The 
metallicity o f hot X-ray gas in clusters has be en measured 
at z ~ 0.3 (|Mu ghj3tszk^_&^JjOc wcnstein 1997 ) and z = 
(Yamashita 1992; Butcher 1995). The metal content of the 
diffuse IGM is estimated from observations of the Lyman-a 
forest (references in figure caption). Figs, [m] and |l^ show 
the model predictions for the average metallicity of the en- 
tire Universe (total mass in metals divided by total mass of 
gas), and for the metallicity of stars, cold gas, hot gas in 
haloes, and diffuse photoionized gas that has been ejected 
from haloes. The cold gas lines from the models may be com- 
pared to the observations of DLAS, the hot gas lines to the 
observations of hot gas in clusters, and the diffuse gas lines 
to the metallicity of the higher-column-density Lyman-a for- 
est. Note that the metals ejected from haloes in our models 
most likely will not escape to large distances, but will re- 
main in the vicinity of the haloes. Therefore our predictions 
for the metallicity of the "diffuse" gas are probably more ap- 
propriately compared with the observational measurements 
in the higher column density Lyman-a forest. 

The differences among the model results are not large. 
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Figure 14. Metallicity as a function of redshift. The square 
symbols are from the measurements o f Zn abundance in DLAS 
([Zn/H DLA ]) from |Pcttini et al. (1997b ). The filled dot at z ~ 0.3 
is fro m measurements of Fc abundance in hot X-ray gas in clus- 
ters ( [Mushotszky fc Loewcnstcin 1997 L and the dot at z = is 



Butcher 



the Fe a bundance from local clusters ( Yamashita 199^ ; ^ 

1995^~^rian gles show estimates of metallicity i n the diffuse IGM 
from ob servations of the Lvman -q fore st (z ~ 0: Shull ct al. 199 



z ~ 0.5: Barlow fc Tvtlcr 1997^ z ~ 3: |R.auch. Haehnelt. fc Stein 
metz 19)7, [goiigaila fc Cowic 1996 , Lu ct al. 1998 , [rytlcr fc Fan 
1994). lines show the collisional starburst model predictions of 



the mass-weighted mean metallicity in the form of stars, cold gas, 
hot gas, or diffuse (photoionized) g labelled on the figure. 



Overall, the collisional starburst model predicts the high- 
est metallicities, the constant efficiency quiescent model the 
lowest. Aside from the DLAS metallicities, which are higher 
than the observations in all three models (see below), the 
largest discrepancy is in the Lyman-a-forest, whose metal- 
licity is too low in the constant efficiency quiescent model. 
The collisional starburst model predicts the flattest metallic- 
ity evolution in all components, implying that the average 
metallicity of stars at z ~ 3 — 4 is already close to solar. 
The very small predicted evolution from z = to z ~ 0.3 
in the metallicity of the hot gas (residing mainly in large, 
cluster-sized haloes) in all three models is in g ood agreement 
with t he observati ons of hot gas in clusters (Mushotszky & 
Loewen stein 1997). 



As noted, the metallicity of cold gas in all models is sys- 
tematically higher than that derived from DLAS. Since the 
metallicities of the cold gas and stars in galaxies are tightly 
coupled, it is difficult to imagine how one could make the 
metallicity of the cold gas much lower yet still match the 
observed solar metallicities of bright galaxies today. On the 
other hand, metallicities estimated from DLA surveys may 
systematically underestimate the true values for several rea- 
sons. First, dusty high-metallicity systems might dim any 
quasar in the line of sight enough t o discourage obs ervers 
from trying to obtain its spectrum (Pei & Fall 1995). Sec- 
ond, the metallicities of the outermost regions of local galax- 
ies are often significantly smaller than their central metal- 



-3 




Figure 15. Same as Fig. [14J for the quiescent model. 



licities. By a simple cross-section argument, there will be 
many more lines of sight passing through these metal-poor 
outskirts than the metal-rich inner regions of the galaxy. Fi- 
nally, in the collisional starburst model, galaxies that have 
experienced a recent burst of star formation, accompanied 
by metal production, will have consumed nearly all of their 
cold gas and will not produce damped systems in absorption. 

An interesting side effect of the collisional starburst sce- 
nario is that it predicts very early, efficient enrichment of the 
IGM and ICM by supernovae driven winds. In our simple 
models, we have no difficulty in ejecting enough metals to 
pollute the IGM to the observed levels by z ~ 3. This is 
consistent with the strong metal-enri ched outflows deduce d 
from the spectra of observed LBGs ( Pettini et al. 1999b ). 
This implies t hat a large fraction of the "missing metals" at 
high redshift ( |Pettini 199s| ; |Pagel 1998|) may be in the form 
of hot gas in clusters or in the IGM. Similarly, enrichment 
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occurs early in the accelerated quiescent model, again be- 



of how sizes are estimated in our models; see Mo, Mao, fa 



cause galaxies with relatively small masses have high star White (199!:) for more detailed modelling of galaxy sizes) 



formation rates. By contrast, in the constant efficiency qui- 
escent model, most star formation takes place in very mas- 
sive haloes, and gas and metals cannot escape from their 
deep potential wells. Significant pollution of the diffuse IGM 
therefore occurs only at much later times (2 < 1). 



Typical values for the mock-catalog galaxies are about a fac- 
tor of two smaller than those of nearby bright galaxies, fur- 
ther evidence that LBGs in the collisional starburst model 
are not fully assembled L*(« = 0) galaxies. The radii of 
the model galaxies are in good agreement with the aver- 
age hal f-light radii of LBGs observed in the H D F (starred 



points; Giavalisco, Steidel, fe Macchetto 1996), 



3.4 PropcrticD of GnlnxicD at 



Lowenthal 



et al. 1997). Model radii at z — are also in good agree- 



Our results so far leave us with the conclusion that the con- 
stant efficiency quiescent model is in serious conflict with 
the observations. The best model seems to be the collisional 
starburst model, although the accelerated quiescent model 
is not strongly ruled out. In this section, we concentrate 
on the collisional starburst model, and investigate the pre- 
dicted properties of z ~ 3 galaxies in detail and compare 
them with the available observations. For this purpose, we 
have created "mock-HDF" catalogs with the same volume as 
the HDF. D ust e xtinction is included using our usual recipe 
(see Section |2.4[ ) . 

Fig. [l^ shows the properties of galaxies in the redshift 
range 2 <C z <C 3.5. The small dots show results as a func- 
tion of the Veoe (rest ~ 1600 A) magnitude, for galaxies in 
a typical mock-HDF catalog. Model galaxies have been se- 
lected with a flat selection function over the redshift range 
2 < 2 < 3.5. The histograms show the distribution of the 
same quantities for galaxies brighter than V606 = 25.5, to 
compare with ground-based samples. Because of the small 
volume of the HDF, we average over many random realiza- 
tions to obtain the histograms. 

Panel a) shows the cold gas fraction / gas = 
cricoid /(n^coid + fistar)- The distribution of / gas is nearly uni- 
form o ver the full range, with a median value of ~ 0.5, and 



ment with the sizes of local galaxies, as we showed in SP. 
The significance of this agreement is less clear here, however, 
since the half-light radii, particularly in the UV, may well be 
considerably smaller than the baryonic half-mass radii that 
we model. However, radii measured in the observed near-IR 
(rest ~4000 A) from the NI CMOS HDF are typica lly quite 
similar to those shown here (Dickinson ct al. 1998). 

Panel d) shows observed and model linewidths. The ve- 
locity dispersions of observed LBGs can be estimated based 
on the widths of emission lines such as H/3 or Om. Emission 
lines have been detected for a few of the brightest LBGs, 
and the velocity dispersions derived from the observed line- 
widths are rather small: 50 — 90 km s -1 for four o bjects, 
and 180 km s" 1 for one object (|Pettini et al. 1998;). The 



widths of the observed emission lines probably reflect the 
mass within about one effective radius. At such small radii, 
there is probably little contribution from dark matter. We 
therefore estimate the 1-D velocity dispersion in our models 
using the expression 



Gm 

cr e 



(•») 



where m is the total baryonic mass (cold gas and stars) 
r e is the effective radius (we use the baryonic half-mass ra- 



dius) , and c is a geometry-dependent factor (Phillips et al 



with many galaxies having values as high as / gas = 0.8 — 0.9. 



This is significantly higher than typical values for local disc 



1997), which we take to equal 6, cor responding to a hot 
comp onent with a density p oc r~ 3 (Binney fa Tremaine 



galaxie; , which arc around 0.1 — 0.25 ( de Blok, McGaugh, 1987). Although there is considerable uncertainty as to how 



fa van der Hulst 1996). We do not have direct observational 
measures of the gas content of the LBG population, but the 
high gas fractions we obtain are consistent with the factor 
of J> 3 increase in the total cold gas content of the Universe 
(^coid) from z — to z ~ 3 from observation s of d amped 
Lyman-Q systems, as demonstrated in Section 3.3.3. 



Panel b) shows the stellar masses of model LBG galax- 
ies. Note that stellar masses range from 10 s to 10 11 Yi~ 2 Mq, 
on average one or two orders of magnitude smaller than the 
stellar masses of present-day L* spirals and ellipticals. Thus, 
according to the collisional starburst model, observed LBGs 
are not the already f ully assembled progenitors of present 



day L > L t galaxies (|Steidel et al. 1996a; Giavalisco, Stei- 



del, fa Macchetto 1996|). It may well be the case that most of 



the stars in these LBGs would end up within bright galaxies 
at z = 0; however, there will be several generations of merg- 
ing in the interim, and no simple one-to-one correspondence 
between the two populations. Accurate stellar masses for 
real LBGs are not yet available; however, the model stellar 
masses that w e obtain are very com patible with the values 
calculated by sawicki & Yee (1998) on the basis of multi- 
band photometry of galaxies in the HDF. 

Panel c) compares the baryonic half-mass radii of model 
galaxies to actual data (see Section 2.3 of SP for a discussion 



to model the linewidths of LBGs, we obtain good agr eement 
with the limited data available ( Pettini et al. 1998| shown 
by the large star symbols). 

In panel a), we remarked on the rather small masses 
of galaxies in the collisional starburst model. That is paral- 
leled here by the modest velocity dispersions. Masses (and 
velocities) tend to be small in this model because small ob- 
jects are elevated to the level of visibility by starbursts. Such 
an effect occurs to a lesser extent in the accelerated quies- 
cent model, but not at all in the constant efficiency quies- 
cent model. This suggests that size indicators of all types 
(masses, radii, velocity dispersions) are a way to discrimi- 
nate among models. In the accelerated quiescent model, for 
example, the shape of the distribution of dispersions is dif- 
ferent — it is flatter and skewed towards larger a v ) but has 
a similar mean — while the CE quiescent model typically 
has the largest dispersions, ranging from 100 to 220 km/s, 
with a mean of about 140 km/s. More detailed modelling 
and additional data will provide important information on 
the baryonic masses of these objects. 

Panel e) presents metallicities for model galaxies. Lit- 
tle is known observationally about observed LBG metallici- 
ties, alth ough wide variations in the strength of Civ absorp- 
tion (see Steidel et al. 1996b, Lowenthal et al. 1997, Trager 
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et al. 1997) suggest that stellar metallicities have a broad 
range of values. The stellar metallicities of our model galax- 
ies are consistent with this, ranging from about one-tenth 
solar to solar, with a mean of about one-half solar. This 
is also quite compatible with the estimate based on inter- 
s tellar absorption line s in the lensed galaxy MS 1512-cB58 
(Pettini et al. 1999b; shown as the dark square on panel 



culiarity and are quite unlikely to resemble classical disc 
galaxies. A quantitative statistical analysis of the morphol- 
ogy of observed LBGs could place important constraints on 
the collisional starburst scenario. We intend to pursue this 
topic in future work (preliminary results are presented in 



Somerville et al. 1999) 



e). A weak metallicity-luminosity relationship is already in 
place at z — 3, with the brightest objects typically showing 
higher metallicities. 

A tight correlation between rest-UV luminosity and in- 
stant aneous star formation rate is often assumed in order to 



Accurate stellar ages of LBGs would provide a good test 
of models. However, the determination of ages from observed 
colours is complicated by the degenerate effect of reddening 
due to dust. Some observational papers discussing the LBG 
population have suggested that th e total duration of star 



formation may b e as large as 1 Gyr ( Steidel et al. 1996b 



estimat a the former from the latter. The solid line in panel 



f) sho ws the relatio n of this sort used by Madau, Pozzetti 



Pet- 



& Dickinson (1997). However, galaxies with complex star 
formation histories, particularly episodic ones, may show a 
non-negligible scatter from this relationship, as shown by 
our model galaxies. Dust extinction also increases the scatter 
in this relation and moves the brightest, most extinguished 
galaxies off of it. The large star symbols in this panel show 



galaxies from the sample of Pettini et al. (1998), where the 



star formation rates have been estimated from ¥L(3 emission 
lines. Star formation rates for the brightest galaxies in the 
collisional starburst model approach 100 solar masses per 
year, in agreement with observed values. The quiescent mod- 
els predict far fewer galaxies with such high star formation 
rates, as reflected in Fig. ^| and ^. 

Panel g) shows the mass ratio of stars in the bulge to 
the total stellar mass (B/T) for the model galaxies. We find 
that bright galaxies are biased towards low B/T ratios, indi- 
cating that the galaxies are disc dominated. Observed B/T 
values (starred points) are also shown for bulge-disc decom- 
positions of HDF galaxie s, using the procedure described in 
Marleau & Simard (1998), where "bulge" and "disc" compo- 



nents were determined by fitting observed surface brightness 
profiles with a Sersic or an exponential form. The availabil- 
ity of photometric redshifts allowed us to select galaxies in 
the redshift range of interest (2.5 < z < 3.5). These obser- 
vations were compiled and provided to us by L. Simard and 
K. Wu. The resulting distribution of B/T with Veoe mag- 
nitude looks very similar to the predictions of our models, 
including the weak apparent correlation of magnitude with 
B/T. However, our B/T values are weighted by mass, while 
observed values are weighted by light. This difference might 
bias the measured values in some unknown way. It is also in- 
teresting to note that about 20 percent of the HDF galaxies 
(brighter than mAs(1500) ~ 26) have B/T > 0.40, similar 
to the fraction of early-type galaxies found in the nearby 
Universe. 

The presence of an exponential profile does not mean 
that a classical disc is present, and similarly an observed r 1//4 
profile does not imply the existence of an early-type galaxy 
in the classical sense. It is clear that traditional morpholog- 
ical definitions must be expanded and modified to usefully 
discuss the inhabitants of the high-redshift Universe. Many 
LBGs show a cen trally concentrated r 1,/4 core within a more 



diffuse envelope ( piavalisco, Stcidel, fc Macchetto 199(: ) , yet 
many also show pronounced substructure and sign s of dis- 



turban ce, even i n the rest- visual NIC MOS images ( Lowen- 
thal et al. 1997; Conselice et al. 1998). The collisional star- 



burst scenario predicts that a measureable fraction of LBGs 
should show significant substructure and morphological pe- 



tini et al. 1997a), based on the TZ — K colours of the objects. 
However, [Sawicki fc Yee (199£ ) conclude that the dominant 



stellar population of the LBGs is less than 0.2 Gyr old, with 
median ages of ~ 10-36 Myr. This conclusion was based on 
their comparison of model SEDs to photometric data from 
5 filter bands (VIJHK) spanning the Balmer break. The IR 
photometric data is helpful in breaking the age-extinction 
degeneracy, although some degeneracy remains. Panel h) in 
Fig. [l^ shows mass-weighted mean stellar ages for model 
galaxies at z = 3 as a function of apparent magnitude (note 
that the luminosity weighted ages, w hich corresponds mor e 
closely to the quantities estimat ed by sawicki & Yee (1998), 
are considerably younger; see § |3.4.1 ). The distribution of 
ages is fairly fiat, with many objects having very young ages 
(indicating that they formed most of their stars in a recent 
burst) and some having ages close to the age of the Universe 
at the mean redshift of the sample rl 



3.4-1 Dust Extinction and Colours 



As described in Section 2.4, we have assumed a very simple 
recipe for dust extinction, in which the optical depth of a 
galaxy is a deterministic function of its intrinsic UV lumi- 
nosity. Recall that the resulting extinction is also a func- 
tion of the galaxy's inclination, which we chose randomly. 
Fig. [It] shows the effects of dust reddening on the model 
galaxies in our mock-HDF catalog as a function of the un- 
extinguished magnitude at rest 1500 A. We convert the ex- 
tinction in magnitudes (shown on the righthand axis) to 
a rest-frame colo ur excess E(B-V) using the recipe given 
in Calzetti (1997). The solid histogram shows the distribu- 
tion for the brightest {attenuated magnitude misoo < 25.5) 
model galaxies. For comparison, we also show the distri- 
bution of rest E(B-V) (provided to us by K. Adelberger) 
obtained from t he spectral slopes o f LBGs from the ground- 
based sam ple of 3teid£l_et_aL_[1999 ) , using the t echnique de- 
scribed by Meurer, Heckman, & Calzetti (199!; ). A detailed 
di scussion of dust extinction i n observed LBGs is presented 



Adelberger & Steidel (2000). Note that negative values of 



E(B-V) are unphysical and result in part from the presence 
of Lyman-a emission in some of the spectra (not included 
in our modelling). In about half of the cases, Lyman-a is 
seen in absorption and this broadens the E(B-V) distribu- 
tion redwards. Within these uncertainties, the distribution 
we obtain in our models is fairly similar to the observed one. 
This gives us some confidence that our simple dust model is 



* The age of the Universe in this cosmology is 3.2 Gyr at z = 2, 
2.1 Gyr at z = 3, and 1.8 Gyr at z = 3.5. 
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Figure 16. Small dots show various properties of galaxies in our mock-HDF catalog, at 2 < z < 3.5 as a function of the Veoe (rest ~ 1500 
A) magnitude of the galaxy. The histograms show the distribution of galaxies in this redshift range that are brighter than Veoe = 25.5. 
The quantities shown are: a) gas fraction, b) stellar mass, c) half-light radius, d) internal velocity dispersion, e) stellar metallicity, f) star 
formation rate, g) bulge-to-total mass ratio, h) mass- weighted mean stellar age. Large star-shaped symbols and the dark square in panel 
e) show observational estimates of some of these quantities (see text). 
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Figure 18. Observed V— I (left) and I— K (right) colours fo r galaxies at z ~ 3. The plus symbols show the colours of galaxies in the 
HDF with photometric redshifts in the range 2.0 < z < 3.5 (jChcn et_al. 1999): JLanzcttact al. 1999ft. Large star symbols show colours 
of HDF galaxies with spectroscopic redshifts ( Lowenthal et al. 1997 ; [gawickT'fc Ycc 199g| ). The small dots show the colours of model 
galaxies extracted from the mock-HDF catalog, assuming a flat selection function in the range 2.0 < z < 3.5. Dust extinction is included 
in the model galaxy colours using our usual recipe. The arrows show the reddening vector for an L* galaxy at z = 3 assuming either a 
Calzetti or SMC extinction curve. 




Figure 19. Two-colour diagram (observed V— I vs. I— K) for galaxies at z ~ 3. In the left panel, large star symbols and plus symbols 
show the colours of HDF galaxies with spectroscopic and photometric redshifts (respectively) in the range 2.0 < z < 3.5 (references in 
text). Small dots show the colours of model galaxies in the same redshift range from the mock-HDF catalog (assuming a flat selection 
function from 2.0 < z < 3.5), including dust extinction as usual. Dot-dashed lines show the unreddened colours of instantaneous burst 
populations of various ages from the solar metallicity Bruzual & Chariot models (GISSEL) at z = 3. The five diamond symbols mark the 
colours of single-burst ages of ~ 25, 50, 100, 250, and 320 Myr, from lower left to upper right. The arrows show the reddening vector for 
an L* galaxy at z = 3, assuming a Calzetti or SMC type extinction curve. The length of the vector represents an extinction of a factor 
of five at 1500 A. 
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Figure 17. Extinction at rest 1500 A for the model galaxies at 
z = 3, as a function of their unattenuated UV magnitude (see 
text for details of the dust recipe). The lefthand axis shows the 
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Figure 20. The (almost) original Madau diagram (Madau 



axis sho vs it in terms of magnitudes of extinction at 1500 A. The 



Pozzetti, & Dickinson 1997). The models shown are: model G 



scatter at fixed magnitude is due to the random inclination chosen 
for each galaxy. The shaded histogram shows the distribution of 
rest E(B-V) for the bright (attenuated misoo < 25.5) galaxies 
in the models. The unshaded histogram shows the distribution 



of rest E(B-V) for LBGs in the ground-b; 



magnitudes by (Adelberger & Steidel 2000) 



,sed sample to similar 



of BCFL, our attempt to reproduce model G of BCFL, and our 
constant efficiency quiescent model. All three are seen to be sim- 
ilar. The dotted line is the curve shown in Fig. 16 of BCFL. Our 
conversion of the data points from luminosity to star formation 
rate is different than the one used by BCFL, and is comparable 
to the star formation in luminous stars, shown for BCFL model 
G by the solid line. 



at least roughly compatible with the best current observa- 
tional estimates. 

We investigate the colours of our model galaxies, includ- 
ing the effects of dust, using the same recipe discussed above. 
Fig. |l8| shows the Veo6 — 7si4 and Isu — K a colours of galax- 
ies in our moc k-HDF (small dots) along with galaxies with 
spectroscopic (Lowenthal et al. 1997; gawicki & Yee 199q) 



and pho tometric redshifts (Chen et al. 1999; Lanzetta et al 



199S). It is immediately apparent that the model galaxies are 
too blup in Ygpg — T slll hy ahnnt. 9 magnitudes while the 



left panel shows a two-colour diagram (Veo6 — -^814 vs. 
J8i4 — K s ) for the observed HDF galaxies with spectroscopic 
(stars) and photometric (crosses) redshifts, as before. It is 
interesting to note the four or five photo-z galaxies with 
anomalously blue 7gi4 — K s colours for their Veoe — hi4 
colours. A cross-check reveals that these galaxies are all 
fainter than 7si4 = 25.5 and therefore are not bri ght enough 
to be contained in the spectroscopic sample of 



7 t. colours are in reasonably good agreement with the et aL [1997 )• By comparing with the dot-dashed lines, which 



Lowenthal 



data. One might be tempted to conclude from this that the 
starburst model produces galaxies with too young a stellar 
population; however, this is unlikely to be the explanation, 
as /8i4 — K a spans the Balmer break at this redshift and 
is a much better indicator of the age of the stellar popula- 
tion. The V606 — 7814 colour reflects the spectral slope in the 
far-UV (~ 1500-2000A), and, as we have discussed, this is 
strongly correlated with the far-IR excess, and hence with 



the amo unt of dust extinction (Meurer, Heckman, & Calzetti 
199S()7]rhis seems to indicate that we have underestimated 
the amount of dust extinction in our model galaxies, or that 
the true extinction curve is steeper than the Calzetti curve 
that we have assumed throughout — perhaps closer to an 
SMC curve. By comparing the reddening vectors shown in 
the figure, one can see that, had we used an SMC curve, 
both sets of colours would have been in quite good agree- 
ment with the data. 

A similar conclusion may be drawn from Fig. [H]. The 



show the unreddened colours of a single-age stellar popula- 
tion from the GISSEL98 models (diamonds mark ages of 
~ 25, 50, 100, 250, and 320 Myr from bottom left to top 
right), we see that their 7gi4 — K a colours are characteristic 
of an extremely young (less than 25 Myr) stellar popula- 
tion but the Veoe — 7si4 colours are considerably redder — 
perhaps these are very heavily extinguished starbursts? The 
righthand panel shows where our model galaxies from the 
mock-HDF lie in this diagram. We do not see these anoma- 
lous objects in the models, even if an SMC extinction law is 
assumed, which could be an indication that our dust mod- 
elling is too simplistic (which would hardly be surprising), or 
that these objects are unusual, have been assigned incorrect 
redshifts, or have large photometric errors. 
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4 COMPARISON WITH PREVIOUS WORK 

The results of this sort of modelling, particularly the red- 
shift evolution, can be quite sensitive to the details of the 
recipes used to model the physical processes. This explains 
some of the differences between our results and those of other 
groups using similar techniques. For example, previous semi- 
analytic modelling of high-redshift galaxies (BCFL) gave a 
rather different picture of the star formation history of the 
Universe and the nature of the Lyman-break galaxies. BCFL 
concluded that observable LBGs form only in very massive 
haloes, > 10 12 /i -1 Mq, and implied that the dominant mode 
of star formation at all redshifts is quiescent. They state ex- 
plicitly, and show in their Fig. 7, that in their models "most 
galaxies... never experience star formation rates in excess of 
a few solar masses per year." This is clearly quite different 
from the results of our best-fitting models. These differences 
arise from the different ingredients that we have chosen to in- 
clude in our models. BCFL made several assumptions that 
made the mass-to-light ratios of their haloes higher than 
ours. Based on our attempts to reproduce their results, we 
believe that some of the important differences are: 

(i) BCFL assumed primordial metallicity for all hot gas, 
reducing the efficiency of gas cooling (see Fig. 2 of SP) . 

(ii) The explicit V c dependence in their star formation 
recipe (Eqn. § of this paper; Eqn. 2.5 and 2.10 of CAFNZ) 
makes star formation less efficient in small V c haloes. This 
suppresses star formation at high redshift, when typical halo 
masses are smaller. 

(iii) The strong Vr dependence of their supernovae feed- 



back r© :ipe (Eqn. 2.8 and 2.11 of CAFNZ) further suppresses 
star formation in small haloes. 

(iv) BCFL assume that only a fraction of the stars are 
luminous. In their fiducial model A, they assume that f* um — 
0.36 (T = [/lunj -1 = 2.8, in their notation), and in model G 
/lum = 0.42. This results in stellar mass-to-light ratios that 
are larger by a similar factor. 

(v) Although merger-driven starbursts are included in the 
models of BCFL, they occur only in major mergers between 
satellites and central galaxies. This neglects the contribution 
from satellite-satellite mergers and minor mergers, which we 
have shown can be quite significant (see Fig. Illl) . 

Fig. shows the star formation history of the Uni- 
verse as represented by the (almost) o riginal Madau diagram 
(Madau, Pozzetti, & Dickinson 1997), along with the results 



of BCFL model G (the closest of their models to our fidu- 
cial cosmology and IMF)|| Also shown are our attempts to 
reproduce their results by adopting the assumptions listed 
above. Many differences in the details of the modelling re- 
main, so it is not surprising that our results do not agree 
exactly, but they are similar enough to give us confidence 
that we understand the main effects. 



t Note that in Fig. 16 of BCFL, the data points shown were 
calculated using a different co nve rsion from luminosity to star 
formatiq 



We also show the results of our constant efficiency quies- 
cent model, the same model shown throughout the previous 
sections. It turns out that the results of model G of BCFL 
are very similar to those of our constant efficiency quiescent 
model. Although the efficiency of quiescent star formation 
is actually lower at early times in the models of BCFL (see 
points (ii) and (iii) above), starbursts in satellite-central ma- 
jor mergers were included (see point (v) above) and this con- 
tributes some additional bright galaxies. One can see from 
Fig. ^ and can further verify by comparing the luminosity 
functions shown in Fig. |s| and Fig. [j] with Fig. 15 of BCFL, 
that, for our purposes, the models produce very similar re- 
sults. Like our CE quiescent model, the models of BCFL 
produced just enough light and just enough bright galaxies 
at z ~ 3 if dust extinction was neglected. We now see that, 
when dust extinction is taken into account, these models 
produce a luminosity function that falls off too steeply on 
the bright end, and predict a steep fall-off in the number den- 
sity of bright galaxies and the integrated luminosity density 
with redshift. This results in very few bright galaxies at red- 
shifts of ;> 5, in sharp contrast to the collisional starburst 
model or the accelerated quiescent model. 

As noted in Section ^, our accelerated quiescent star 
formation recipe is effectively very similar to the one usu- 
ally used by the Munich group (e.g. KWG93). We have 
shown that this recipe produces fairly good agreement with 
many of the observations, but does not produce enough 
very bright LBGs and also underpredicts the amount of 
cold gas at high redshift compared with o bservations of 
damped Lyman-a systems. In a recent paper, Kauffmann & 



Haehnelt (2000) reached a similar conclusion about the cold 



gas problem and adopted a similar solution. They write their 
quiescent star formation recipe as rh* = am co id/tdyn- As 
noted before, the dynamical time td yn scales approximately 
as tdyn oc (1 + z)~ z ^ 2 , so their adopted scaling of the param- 
eter a with redshift as a(z) = q(0)(1 + z) -7 with 7 = 1-2 
is very similar to our "constant efficiency quiescent" recipe. 
They also include bursts in major satellite-central mergers, 
and find that these bursts are responsible for about two- 
thirds of the star formation at high redshift. They find that 
this model, which is similar to our favored collisional star- 
burst model, produces good results for the redshift evolution 
of the space density of bright quasars. Taken together, this 
suggests a picture in which the brightest inhabitants of the 
high-redshift Universe, LBGs and quasars, may be produced 
by the same process: strong inflows induced by mergers. 

It is also interesting to compare our results to those 
of cosmological A^-bod y simulations including hydrodynam - 
ics and star formation. [Katz, Hernquist, fc Weinberg (19"9£ ) 
performed a detailed study of the clustering of high red- 
shift galaxies in such simulations, and concluded that the 
observed clustering and number densities could be easily 
explained in a variety of CDM models. In this respect, their 
results are quite similar to the massive halo models discussed 
in Section |l| However, although their simulations included 
gas cooling and star formation, they assumed a monotonic 
relation between baryonic mass and UV luminosity. As we 



n than oura (3cc Tabic in our Appendix, and Tabic 5 

Of BCFlj la addition, RflFT, plotted "tntal" atav f/",rmafir,n 



have already discussed, a nd as acknowledged by Katz, Hern 



quist, fe Weinberg (1999 ) , the effects of episodic star forma- 



rate in the models and corrected the observations assuming that 
only a fraction [T] _1 of the stars are luminous. We instead plot 
the star formation of luminous stars in all of the models. 



tion and dust could significantly complicate this relation- 
ship. A more direct look at the na ture of the z ~ 3 g alaxies 
in these simulations was taken by Dave et al. (199!: ) . Here, 
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the star formation history of the galaxies was extracted and 
convolved with stellar population synthesis models to ob- 
tain estimates of their luminosities, and the effect of dust 
was considered. They then find that the luminosity function 
at 2 ~ 3 is in agreement with the observations when dust ex- 
tinction is included. The galaxies brighter than 25.5 in their 
simulations have fairly large stellar masses (log m sta r > 10.4 
for their value of /iioo=0.65) and contain a significant older 
stellar population; i.e. they are not predominantly young 
bursts. However, Dave et al. acknowledge that their simula- 
tions do not have sufficient resolution to properly model the 
collisional starburst mechanism that we have invoked here, 
so they cannot rule out the possibility that small starbursts 
could also be present. It is also worth noting that in these 
simulations, the star formation efficiency scales as the gas 
density to some power, which will give a similar behaviour 
to our "accelerated quiescent" star formation recipe, and 
this may partially explain why they obtain sufficient num- 
bers of LBGs without the need for an additional population 
of starbursts. The volume of their simulations is extremely 
small (~ ll/i -1 Mpc on a side), far too small to probe the 
very bright end of the luminosity function, where we found 
the most pronounced differences between the collisional star- 
burst model and the accelerated quiescent model. 



5 DEPENDENCE ON MODELLING 
ASSUMPTIONS 

This se ction discusses the sensitivity of our results to the 
cosmology, IMF, and stellar population modelling. 
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Figure 21. The star formation history in an Q = 1 universe. 
Small open and large filled points show the observations with- 
out and with (respectively) corrections for dust extinction (see 
Appendix). The solid line shows a collisional starburst model in 
a tCDM cosmology, normalized to the Tully-Fisher relation at 
2 = 0. The dotted line shows the same model norm alized to the 
prese nt-day dust-corrected UV luminosity density (Treyer et al 



1998) 



5.1 Cosmology and Power Spectrum 

With improving observational constraints, the freedom to 
choose cosmological parameters as one wishes is rapidly dis- 
sappearing, and therefore we have chosen to focus on the 
popular ACDM cosmology. Any cosmology within the cur- 
rent observationally favoured regime of parameter space, 
namely £1 — 0.3 — 0.5 (open or flat) and F ~ 0.2, would 
give similar results. The old standard, cluster normalized 
SCDM, would also give similar results because of the large 
amount of power on small scales. Models with = 1 and 
realistic power spectra, such as tCDM, CHDM, or tCDM 
(see SP) show a stronger peak at z ~ 1.5 and a sharp de- 
cline in the star formation rate at higher redshift. In Fig. ^ 
we show the star formation history for a model in which we 
use the rCDM cosmology defined in SP and normalize to 
the local Tully-Fisher relation in our usual way. This figure 
should be compared to the corresponding Fig. ^| and [| for 
our standard ACDM cosmology. The rCDM model overpro- 
duces the luminosity density at low redshift (z <J 2) and 
falls below the observations at high redshift (z ~ 4). If we 
normalize the model to fit the local UV-luminosity density, 
we find that there is not enough star formation at high red- 
shift compared to the dust-corrected new Madau diagram. 
These models also violate the constraint on O gas from DLAS 
at high redshift. The results for models with similar power 
spectra such as CHDM or tCDM would be nearly identical. 
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Figure 22. The luminosity of a single burst of mass 10 6 Mq as 
a function of age for the GISSEL models with solar metallicity 
and various IMF. Retaining the Salpeter shape but changing the 
lower mass cutoff from 0.1 Mq to 1 Mq would result in shifting 
the luminosities up by about a factor of three. 
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Figure 23. Same as Fig. |22|, but for the STARBURST99 mod- 
els with different metallicities. The far-UV luminosity does not 
depend much on the metallicity for the very young stars that 
dominate it. 



5.2 Stellar Population Synthesis 

Most of the currently available observations of galaxies at 
z ;> 3 are viewing light that is emitted in the far UV (~ 
1500 A). This light is dominated by very young (< 10 Myr), 
massive O and B type stars, which are notoriously poorly 
understood and difficult to model. In addition, the stellar 
mass-to-light ratio depends on the metallicity and IMF. 

We have made a comparison of the luminosity at 1500 
A (averaged over a tophat with width 400 A) of an instanta- 
neous burst for a variety of the stellar populatio n synthesis 
models that are publicall y available: GISSEL (Bruzual & 
Charlol | 199^ ) PEGASE ([Fioc fc Rocca-Volmerange 1997|) , 
STARDUST flDcvriendt, Guiderdoni, fe Sadat 1999[ ), and 

STARBURST99 ( [Leitherer et al. 1999| ) p|. For models with 
solar metallicity and a Salpeter IMF, all of these models 
agree within 0.15 dex for populations younger than 10 Myr 
old, which dominate the UV light especially in starburst 
galaxies. The GISSEL, STARDUST, and PEGASE models 
agree at a similar level for older populations as well, but the 
STARBURST99 models differ from the others by as much 
as a factor of fifteen for populations up to 1 Gyr in age. We 
find a similar level of agreement for populations with 0.2 
and 0.4 solar metallicity. 

Fig. ^ shows the dependence on the assumed IMF for 
the solar metallicity GISSEL models. We see that the UV 
luminosity is quite sensitive to the IMF, in particular to 



t The GISSEL models were provided to us by S. Chariot, the 
PEGASE models by M. Fioc, and the STARDUST models by J. 
Devriendt. The S TARBURST99 models were downloaded fro m 
the public website http://www.stsci.edu/scicnce/starburst99/ 



the fraction of high-mass vs. low-mass stars. The Scalo IMF 
gives results that are lower by about a factor of three at 
ages younger than 10 Myr. Similarly, changing the lower- or 
upper-mass cutoff would have an effect on the results. For 
example, if we drastically increased the lower-mass cutoff to 
1 Mq instead of 0.1 Mq and retained the Salpeter shape, the 
curve would shift up by about a factor of three |§J None of 
these changes are huge in the context of the theory, but evo- 
lutionary changes in the IMF as a function of redshift could 
induce a tilt in the Madau diagram, and thus affect com- 
parison with observations. However, to salvage the constant 
efficiency quiescent model would require a total evolutionary 
change of a factor of 10 in mass-to-light ratio (see Fig. [|), 
which would stretch all the above effects to their limits. 

The STARBURST99 models have been constructed to 
be particularly relevant to the very young stellar popula- 
tions that we are interested in here, and use different stellar 
atmosphere models and a different approach to modelling 
mass loss than the ot her three models mentioned above (see 
Leitherer et al. 1999[ ). In Fig. ||, we show the dependence 
on metallicity in these models. For models with metallicities 
ranging from one-fifth solar to solar metallicity, the lumi- 
nosity at 1500 A is almost independent of metallicity for 
stars younger than 10 Myr. The luminosity becomes more 
sensitive to metallicity for older stars (;> 200 Myr). 

Note that we have neglected the contribution from neb- 
ular emission in our modelling. In the STARBURST99 mod- 
els (Salpeter with solar metallicity) nebular emission con- 
tributes about 15 percent to the luminosity at 1500 A at an 
age of 1 Myr. 



5.3 Dust and Reddening 

Our treatment of dust and reddening is perhaps the least 
secure of our assumptions. This treatment falls into three 
broad categories: the mean absorption per typical galaxy, 
the variation in absorption versus wavelength (reddening 
curve), and the variation in absorption versus galaxy lu- 
minosity (brighter galaxies are more extinguished). None of 
these threaten the two main conclusions of this paper, that 
the collisional starburst model is consistent with present 
data on LBGs, and that star formation efficiency at high 
redshift must be higher than locally. With regard to the 
first conclusion, we are claiming only consistency using the 
best information presently available. This has certainly been 
shown. With regard to the second point, even if we reduce 
the the dust correction to zero, the constant efficiency qui- 



escent model is inconsistent with the Steidel et al. (199!:) 
data at z ~ 4 and preliminary detections of z :> 5 galaxies. 

There are hints that some aspects of our dust treatment 
are in fact in error. For example, our model galaxies are 
systematically 0.2 mag too blue in V— I (Fig. [l8]). Fixing this 
may necessitate changing the reddening curve to a steeper 
one closer in shape to the SMC curve; however, the overall 
amount of absorption would still remain roughly the same 
(see figure). Our third assumption (that brighter galaxies are 
more extinguished) is based in part on precisely the HDF 
galaxies we are modelling — after correction for extinction, 



3 We base this on a model that we ran using the STARBURST99 
website 
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is it preci sely the reddest objects that are b olometrically 
brighte st ( Meurer, Heckman. fc Calzctti 1999|). 



Opinion seems to be converging on our first and 



most crucial assumption, the average absorption per typical 
galaxy. Given the observed V— I colors of LBG galaxies, it 
is inconceivable that they have zero reddening. At the other 
extreme, dust absorptions three times larger than we have 
assumed w ould predict many more sub-millim eter sources 
than seen (Meurer, Heckman, & Calzetti 1999). Tweaking 



the mean absorption within these boundaries would not 
greatly alter our picture. 

On the other hand, a serious change to our assumed 
absorption treatment would be the discovery that mean ab- 
sorption at a fixed luminosity declines strongly with redshift 
(the reverse variation is implausible). This would introduce 
a strong bias into all tests versus time, which is the most im- 
portant probe of LBG models. Present data on the very dis- 
tant Universe are simply too sparse to know for sure whether 
distant sources are systematically less obscured. Future sub- 
millimeter observations are needed to pin down this last, 
very important point. 



6 SU MMARY AND CONCLUSIONS 



The process of star formation is one of the largest uncertain- 
ties in galaxy formation modelling. We have explored three 
different recipes for star formation and the implications for 
observations of galaxies at high redshift. We investigated 
a model in which quiescent star formation has a constant 
efficiency, and galaxy-galaxy mergers trigger bursts of star 
formation. In this model, we find that most of the observ- 
able Lyman-break galaxies are collisional starbursts. We also 
considered two models in which star formation occurs only 
in a quiescent mode, either with constant or "accelerated" 
efficiency. 

We find that our collisional starburst model produces 
excellent agreement with the observed number density of 
bright galaxies as a function of redshift from 2 < z < 6 when 
a reasonable amount of dust extinction is accounted for. The 
accelerated quiescent model produces nearly the same num- 
ber density of bright galaxies as the starburst model over 
this re dshift range. In contrast to the constant efficiency 



fects of dust into the models using a simple, e mpirically 
motivated mod el of differential dust extinction (Wang & 



Heckman 1996), we obtain very good agreement with the 



actual observed luminosity function. Both the accelerated 
quiescent and the constant efficiency quiescent models show 
a steeper decline at the bright end. The constant efficiency 
quiescent model produces reasonable agreement with the ob- 
served bright end of the luminosity function at z ~ 3 if dust 
extinction is neglected, but with a realistic correction for 
dust extinction, these models predict a cutoff on the bright 
end that is much too steep. The accelerated quiescent model 
also has the wrong sh ape compared to the dust-co rrected 
luminosity function of Adelberger & Steidel (200C). This 
model therefore underpredicts the number density of the 
very brightest LBGs (1Z ab < 22 at z ~ 3). 

The global star formation history of the Universe is re- 
flected in the redshift evolution of the total densities of star 
formation, cold gas, stars, and metals. We have compiled 
a new Madau diagram, using recent observations and an 
observationally based dust correction. We find a broad con- 
sistency between the dust-corrected optical estimates of the 
star formation rate density and tho se obtained from f ar-IR 
observations at z ~ 0.7 from ISO (Flores et al. 1999) and 
sub- mm observations at z ~ 3 from SCUBA (Hughes et al 



1998j ). The shape of the new Madau diagram is quite differ- 



ent from the original one, with a more gradual rise in the star 
formation rate from z = to z ~ 2, and a plateau thereafter. 
The shape of the Madau diagram produced by our models is 
quite sensitive to our adopted recipe for star formation. The 
collisional starburst model is in excellent agreement with 
the new dust-corrected data, and the accelerated quiescent 
model is also in reasonable agreement. The constant effi- 
ciency quiescent model shows a much steeper decline with 
redshift, and have too little star formation at high redshift 
compared with the new data. 

Similarly, the collisional starburst model and the ac- 
celerated quiescent model form their stars at a fairly high 
redshift, and thus stellar ages are in agreement with obser- 
vational constraints from "fossi l evidence" co ntained in the 
ages of stars in nearby galaxies ( Renzini 1998] ) . The constant 
efficiency quiescent model forms its stars too late. In addi- 
tion, we show that the redshift evolution of the de nsity of 



cold gas is consistent with constrain ts from DLAS (Btorric- 



quiesce jit model, which predict that very few bright galax- Lombardi, McMahon, fc Irwin 1996|) in the collisional star- 



ies would be in place at z > 4, the burst model and the 
accelerated quiescent model make the tantalizing prediction 
that the comoving number density of bright galaxies at red- 
shift 5 and 6 is nearly as high as at z ~ 3. While this is 
consistent with the numbers of galaxies in the HDF-N and 



S with photometric r edshifts in this range (Chen et al. 1999 



Lanzetta et al. 1999), and some spectroscopic confirmations 



of very high redshift galaxies exist, secure estimates of the 
number densities of galaxies at these very high redshifts may 
have to wait for NGST and results from the new generation 
of infrared spectrographs on large ground-based telescopes. 

A characteristic feature of the collisional starburst 
model is that the unextinguished rest-UV luminosity func- 
tion at z ~ 3 is much flatter than a Schechter function at 
the bri ght end. Our prediction of the intrinsic luminosity 



burst model. The constant efficiency quiescent model may 
have too much cold gas at high redshift, while the acceler- 
ated quiescent model may have too little. 

We investigate the predicted metal enrichment history 
of the Universe, and find that an interesting side-effect of the 
burst and accelerated quiescent scenarios is efficient pollu- 
tion of the IGM and ICM at high redshift due to supernovae- 
driven outflows. This may help to explain the large cover- 
ing factor implied by the high probability with which Civ 
and Mgn absorption are seen in QSO spectra whenever a 
bright galaxy at the ap propriate red shift is < 30/i _1 kpc of 
the QSO line of sight ( |Steidel 1995 ). It also implies that a 
substantial fraction of the metals at high redshift (z ~ 3) 
are in the form of hot gas in haloes or in the diffuse I GM, 



suggestin g a solution t o the "missing metals" problem (Pet- 



functio: l is in excellent agreement with the dust-corrected 



tini 1999 



luminosity function of observed LBGs calculated by Adel- 
berger fc Steidel (200C). Conversely, if we introduce the ef- 



Pagel 1998). When normalized to produce solar 



metallicity stars in bright galaxies at redshift zero, with no 
introduction of additional free parameters, the burst model 
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simultaneously predicts the correct observed metallicity of 
diffuse gas (Lyman-a forest) at z ~ 3, as well as the observed 
metallicity of hot gas in clusters at z = — 0.3. It is natural 
to associate the cold interstellar gas in our model galaxies 
with the observed population of DLAS. Here an interesting 
discrepancy is that the cold gas in our models has a system- 
atically higher metallicity than the DLAS, by about a factor 
of ~ 3 at all redshifts. Moreover, as the constraints on the 
metallicities of very low redshift (z < 1.0) DLAS improve 
( Pettini et al. 1999a ), if this gas is indeed the reservoir for 
star formation, it becomes difficult to understand how it can 
become enriched to solar values over such a short time inter- 
val. This favours the explanation that the observed popula- 
tion of DLAS systematically underestimates the metallicity 
of the total cold gas in galaxies. 

The collisional starburst model produces good agree- 
ment with observable individual properties of z ~ 3 galaxies, 
such as the half-light radii and internal velocity dispersions. 
Our models also predict that high redshift galaxies should 
have cold gas fractions that are much larger than present day 
galaxies: values of / gas = m sas /(m t +m gas ) of 0.5 to 0.9 are 
typical. The model galaxies at z ~ 3 have perhaps rather 
surprisingly high metallicities, from 1/3 solar to solar for 
bright galaxies. The stellar masses of the model LBGs range 
from 10 s — 10 11 h -2 Mq, up to three orders of magnitude 
smaller than present-day L* galaxies. If our models reflect 
the real Universe, this indicates that the majority of the 
LBGs are not the fully-formed, direct progenitors of today's 
L» galaxies, but must experience considerable growth (via 
accretion or merging) by z — if they are indeed the pro- 
genitors of massive present-day galaxies. Another possible 
destiny for some of the smaller-mass LBGs 10 9 h" 2 M Q ) 
is that their stars, gas, and globular clusters are stripped 
as they fall into the potential well of a nearby massive dark 
matter halo, for ming a Pop II ste llar halo such as that of 
the Milky Way (Trager et al. 1997). This resembles the pic- 
ture of g alactic stellar halo for mation outlined in the classic 
paper by Searle fc Zinn (1978), which is supported by much 
recent evidence (see Majcwski, Munn, fe Ha wley 1996, Car- 



ney et il. 199C, Bommer-Larsen et al. 1997, and references 
therein). 

Finally, we compared our results with previous work 
and discussed the effects that different model assumptions 
have on our results. In the models of BCFL, star forma- 
tion was made inefficient in small galaxies by the combined 
properties of the recipes for star formation and supernovae 
feedback. This led to the suppression of star formation at 
high redshift, and a prediction of a strongly decreasing star 
formation rate density with redshift. We argued that our 
constant efficiency quiescent model gives overall very sim- 
ilar results to the BCFL models, and therefore the same 
conclusions will apply: namely, these models are not con- 
sistent with the recent data at z ~ 3 and z ~ 4 when 
the observationally favored correction for dust extinction is 
taken into account. The star formation recipe usually used 
by Kauffmann et al. is very similar to our a ccelerated qui- 
escent model. It is interesting however that Kauffmann & 



Haehndt (2000) recently exchanged this recipe in favor of 



one that is effectively rather similar to our model with con- 
stant efficiency quiescent star formation plus collisional star- 
bursts. They disfavored the "accelerated" type model both 
because they found, as we did, that too much gas was con- 



sumed to be consistant with the observations of DLAS at 
z J> 2, but also because it did not reproduce the observed 
redshift evolution of the space density of bright quasars. 
Thus the same mechanism, merger-driven inflows, may ac- 
count for both high redshift galaxies and quasars. 

We showed that in cosmologies with = 1 and realistic 
power spectra (e.g., rCDM, CHDM, tilted CDM), the de- 
cline in the star formation rate at high redshift is too steep 
to be consistant with the dust-corrected data even in the 
burst model. However, in any cosmology with parameters 
close to the values favored by a broad range of observations 
(Oo = 0.3 — 0.5), we would have obtained similar results 
to those presented here. We showed a comparison of the 
UV stellar-mass-to-light ratio in a variety of different stel- 
lar population models, and showed that there is reasonably 
good agreement between the models produced by different 
groups. We showed that the UV luminosites would have been 
about a factor of three lower had we assumed a Scalo instead 
of Salpeter IMF, and about a factor of three higher had we 
assumed a larger value for the lower mass cutoff (1 Mq in- 
stead of 0.1 Mq). We also showed that according to these 
models, the UV luminosity in galaxies with active star for- 
mation is nearly independent of the metallicity of the stellar 
population. 

We summarize our main conclusions as follows: 

• The details of the recipes used to model star formation 
can have very large effects on the results of galaxy formation 
models (either semi-analytic or numerical), especially when 
redshift evolution is considered. 

• Models in which the efficiency of star formation is con- 
stant with redshift are strongly inconsistent with observa- 
tions at high redshift (z > 3) when the effects of dust ex- 
tinction are taken into account, and are inconsistent with 
higher redshift observations (z :> 4) even when dust is ne- 
glected. Thus, a rather robust conclusion of this study is 
that the efficiency of star formation (i.e. the star formation 
rate per unit mass of cold gas in a galaxy) must increase 
with redshift. 

• This increased efficiency can be accomplished in at least 
two physically plausible ways; either due to collisional star- 
bursts or to the scaling of the star formation rate with dy- 
namical time or gas surface density. Both models produce 
good agreement with most of the observations considered 
here. 

• We favor the collisional starburst mechanism because 
it gives better agreement with the shape of the z ~ 3 lu- 
minosity function at the very bright end, and it is in better 
agreement with the constraints on the density of cold gas 
at high redshift from DLAS. However, due to the remain- 
ing uncertainties in the modelling, we do not consider the 
accelerated quiescent model to be strongly ruled out. 

The two successful models are based on completely dif- 
ferent physical ideas but are difficult to distinguish conclu- 
sively from the present observations. It is crucial to eventu- 
ally determine which process is actually the dominant mode 
of star formation at various redshifts. In one case, by study- 
ing high redshift galaxies we can expect to learn something 
about the merger rate and the efficiency of merger-driven 
inflows. In the other case, we expect to learn more about 
the masses and internal properties of discs. 

Both scenarios have observational support and theoret- 
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ical motivation. It may of course be that the true situation 
involves a combination of both scenarios, however, simply 
adding starbursts to our usual accelerated quiescent recipe 
does not provide a solution. This is because the accelerated 



Table Al. Conversion factors from luminosity to star formation 
(SFR = C X L), where SFR is the star formation rate in MQyr -1 
and L is the luminosity in erg/s for Ha. and the luminosity density 



quiesce nt star formation very rapidly consumes the cold gas 
supply at high redshift. The contribution from the burst 
mode remains small because of the decreased gas fractions, 
and the constraints on the cold gas at high redshift from 
DLAS are badly violated. Real progress in this sort of mod- 
elling will only be made by replacing some of our very simple 
recipes with more detailed and physically motivated pre- 
scriptions. This will require a better understanding of star 
formation in both normal and starburst galaxies. 

An important observational aspect of the LBGs, which 
we have not addressed at all in this paper, is their clustering 
properties. Clustering in the collisional starburst scenario 
cannot be accurately modelled with analytic approaches. 
Therefore we study this using high-resolution N -body simu- 
lations in a companion paper ( ]Koratt et al. 1999| ). One might 
hope to discriminate between the collisional starburst and 
accelerated quiescent scenarios using the clustering proper- 
ties of galaxies, particularly the numbe r of close pairs. How - 
ever, we show in a forthcoming paper (Wechsler et al. 200C) 
that this is not possible. We find that both scenarios are 
consistent with the observed clustering properties at z ~ 3. 

There are, however, some direct observational tests 
which are feasible in the near future and which may begin to 
discriminate between the two scenarios. The collisional star- 
burst scenario predicts that there should be a relatively large 
population of bright, heavily extinguished galaxies with star 
formation rates of hundreds of solar masses per year. There 
are already preliminary indications from the SCUBA results 
that this population has been detected. Future sub-mm ex- 
periments with higher resolution and sensitivity will put 
stronger constraints on the actual numbers and the asso- 
ciated total star formation rates associated with these ob- 
jects. Another way of distinguishing the scenarios is from the 
morphologies of galaxies at high redshift. In the collisional 
starburst picture, we expect many of the LBGs to appear 
highly disturbed and to contain significant substructure. In 
the accelerated quiescent scenario, we expect the galaxies to 
appear smaller and denser, but otherwise similar to normal 
local spirals. Again, there are preliminary indications that a 
large fraction of the observed LBGs do show str ong distor- 



tion an d sub-structure, even in re st visual bands (Dickinson 
Conselice et al. 1998), but a more quantitative 



et al. 1 )98 



analysis is needed. Finally, the observational sample of LBGs 
with measured emission linewidths is growing rapidly. With 
better statistics and more detailed modelling, these data will 
also help to discriminate between the two scenarios. 
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APPENDIX A: HOW TO DRAW A MADAU 
DIAGRAM 

The star formation rate per unit comoving volume as a func- 
tion of redshift was first compiled from observa tions extend- 
ing fro m redshift zero to a redshift of about 4 by 



Vladau et al 

( 1996j ) . The now- famous "Madau diagram" sketched out a 
picture of the history of star formation, from a very early 
epoch when galaxies were perhaps first forming until the 
present day. It has become popular to add more and more 
points to this diagram. However, unfortunately, as different 
authors have added their own points, the calculation of the 
derived quantity, the total star formation rate density, from 
what is actually observed (luminosities of galaxies selected 
in some way) has not always been consistent. We therefore 
think it is timely to revisit the steps in calculating this de- 
rived quantity from the observations, and to compile a set 
of points that have been calculated in as consistent a way 
as possible. As most of the results quoted in the literature 
assume an Einstein de Sitter cosmology, we also provide the 
conversion to other cosmologies. 

To draw your own Madau plot, follow these steps: 

(i) Correct for Incompleteness 

Observational samples are generally flux limited, and thus 
the intrinsic luminosity of the faintest objects in the sample 
changes with redshift. In order to understand the true red- 
shift dependence of the total luminosity density, one must 
correct the observations for incompleteness. This is most eas- 
ily done by fitting a functional form (i.e. a Schechter func- 
tion) to the luminosity function obtained from the observa- 
tions themselves. If the usual parameters of the Schechter 
function, (/>*, L*, and a, are given, then the total luminosity 
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Table A2. Incompleteness corrected luminosity densities at various redshifts, using various tracers of star formation, for Q = 1 (go = 0.5; 
EDS), Co = 0.3, Q A = 0.7 (ACDM) and H = 0.3, C A = 0.0 (OCDM). Units are h" 2 erg s" 1 Hz" 1 . 



Reference 



L 1 clLv. 1 


V 1 1 f 1 Q 1 1 1 f +~ 


1 r EDS] 




1 r OCDMi 


1 OCT f^WClV 


Ha 
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-1.80 


-1.80 


-1.80 


0.2 


Ha 


0.0425 


-1.60 


-1.62 


-1.61 


0.05 


£2000 


0.15 


-1.70 


-1.76 


-1.72 


0.02 


Ha 


0.2 


-1.46 


-1.53 


-1.49 


0.04 


^2800 


0.35 


-1.71 


-1.83 


-1.75 


0.07 




0.625 


-1.39 


-1.56 


-1.46 


0.08 
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-1.15 


0.15 


^FIR 
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-0.93 


-1.10 


-0.98 


0.08 


£<2000 


0.7 


-1.31 


-1.49 


-1.39 


0.06 




1.25 


-1.11 


-1.33 


-1.22 


0.15 


^2800 


0.75 


-1.08 


-1.26 


-1.16 


0.15 
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-0.91 


-1.13 


-1.02 


0.15 
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-1.14 


0.15 


^1500 
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-1.4 
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0.15 
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0.15 
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0.2 
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0.07 
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0.16 


^2800 
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0.04 
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-1.40 
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(+0.31)(-0.22) 




0.75 


-1.17 


-1.35 


-1.25 


(+0.23)(-0.14) 




1.25 


-0.98 


-1.20 


-1.09 


(+0.24)(-0.12) 




1.75 


-0.95 


-1.19 


-1.08 


(+0.24)(-0.12) 




2.5 


-1.16 


-1.41 


-1.31 


(+0.28) (-0.21) 




3.5 


-1.19 


-1.44 


-1.36 


(+0.34) (-0.27) 




4.5 


-1.08 


-1.34 


-1.26 


(+0.44) (-0.37) 
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-1.44 
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-1.63 


(+0.59)(-0.38) 
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Flores et al. 1999 
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density is given by (j> t L,r(2 + a). Here F is the usual Gamma 
function (Beyer 1987). Note that for values of a steeper than 
— 1, faint galaxies contribute a substantial fraction of the 
total luminosity density and therefore the results are quite 
sensitive to the faint end slope, which is often poorly con- 
strained. This step leads to the first source of inconsistency: 
different authors have assumed different values of a, or have 
integrated down to different lower limiting luminosities, 
(ii) Convert to the Desired Cosmology 
Most observational references quote luminosity densities 
assuming an Einstein de Sitter (Q — 1) cosmology. To con- 
vert from one cosmology to another, one must take into 
account two effects. The luminosity derived from a given 
apparent magnitude will change, as will the comoving vol- 
ume derived from a given angular size and redshift range. 
Luminosities scale as 



log ( \ 



2 log 



where L 



old/n 



df\z) 
dT w (z) 



and ef^ ld ^ now (z) are the luminosity and lu- 



(Al) 



minosity distance at a given redshift in the old and new 
cosmologies. If we now define fv = V n ew/V \d as the ratio of 
the comoving volume in the new cosmology to that in the 
old cosmology, then the luminosity density scales as 



log 



2 log 



,/ncw 



■logLM*)]. 



(A2) 



(iii) Convert Luminosity to Star Formation 

The usual tracers of star formation are the luminosity of 
nebular emission lines like Ha or On, or the far-UV (1500- 
2800 A) continuum. This observable quantity must be con- 
verted into a star formation rate. This conversion generally 
relies on stellar population models an d an assumed star for- 
mation history and IMF (see Section 5.2 for a discussion of 



the attendant uncertainties and a quantitative comparison 
of different stellar population models). We give a compila- 
tion of conversion factors for various tracers of star forma- 



tion in Table Al. These are taken from Madau, Pozzetti, & 
Dickinson (1997| ), assuming a Salpeter IMF. The use of dif- 
ferent conversion factors is a second source of inconsistency 
in published results in the literature. 

(iv) Correct for Dust Extinction 

If the tracer of star formation is an optical or UV luminos- 
ity, then the effects of dust extinction may be non-negligible. 
In the original Madau plot, no attempt was made to correct 
for extinction due to the lack of knowledge at that time 
about the effect that dust was likely to have on the obser- 
vations, particularly at high redshift. As we discussed in the 



main text (see Section 2.4), some observational estimates 
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of the extinction in the UV (~ 1500-2000 A) are now avail- 
able. These seem to indicate that the amount of extinction 
at this wavelength in nearby starburst galaxies is similar to 
that in the z ~ 3 LBGs, if the correlation between spectral 
slope and dust extinction remains the same. However, the 
amount of extinction depends strongly on luminosity, with 
instrinsically brighter (more rapidly star-forming) galaxies 
being more heavily extinguished. Therefore it is probably 
not a good idea to apply a fixed correction factor to galax- 
ies of all luminosities, as we would in effect be doing if we 
first corrected the luminosity density for incompleteness, as 
described in step 1, and then multiplied by a fixed factor 
to correct for dust extinction. If this correction factor was 



Steidel et al. 1999). Note that according to our previous cal- 



derived from observations of bright galaxies, but a substan- 
tial fraction of the total luminosity density is contributed by 
faint galaxies, this will lead to an overestimate of the total, 
dust-corrected luminosity density. 

We give a compilation of star formation densities de- 



rived from various tracers at various redshifts in Table A2 
for three different cosmologies (see Table caption). These 
values have been corrected for incompleteness by integrat- 
ing over the entire luminosity function and converted from 
lumino sity to star formation rate using the conversion fac- 



tors fr qm Table Al , but no correction for dust extinction 
has been made. Fig. ^ shows the resulting Madau plot for 
the ACDM cosmology, and Fig. ^ for an Slo = 1 cosmology. 

We now select a subset of these observations, those 
which we believe to give the most robust estimates of the 
star formation rate density over a broad range of redshifts. 
These samples represent a reasonably large volume (unlike 



the res ults from the HDF), are based on spectroscopic red- 
shifts (unlike the results based on the more uncertain pho- 
tometric redshifts), and were observed in roughly the same 
rest waveband (1500-2000 A), and so do not require large 
photometric extrapolations or interpolations. We calculate 
the integrated UV luminosity den sity at l ow redshift from 
the results of Treyer et al. (1998), from Oowie, Songaila 
& Bar jer (1999^ " at intermediate redshift, and from [gtei- 



del et a l. (1999[ ) at high redshift. The sample with the best 
constrained faint-en d slope is the lowest redshift sample of 
Treyer et al. (1998), which has a derived slope a — —1.6, 
exactly the same as the faint end slope at z = 3 derived 
f rom the combined ground-based and HDF samples of LBGs 
( steidel et al. 1999 ) . We see no particular reason that the U V 



luminosity function should flatten at intermediate redshifts, 
so we assume a = —1.6 for the Cowie et al. sample as well, 
and calculate all the incompleteness correc tions accordingly 
(we us e the valu es of 0, and L» quoted by Cowie, Songaila, 
& Barg jer (1999]) for a fixed a = —1.5, but we actually use 
ct = — 1.6 in computing the luminosity density). 

The results of Steidel et al. (1999) indicate that there 



is a factor of 4.7 extinction at ~ 1500 A in LBGs brighter 
than 1Z = 25.5, which corresponds to 0.4L, at the mean 
redshift of the sample (z ~ 3) and for our assumed cosmol- 
ogy. We will assume that these results hold for all of the 
UV-selected samples at all redshifts. To obtain a conserva- 
tive ("minimal") dust correction, we apply the factor of 4.7 
correction to all galaxies brighter than 0.4L, (using the ap- 
propriate value of L* at each redshift) and no correction for 
galaxies fainter that this limit, for each of the three sam ples 
( Treyer et al. 1998 , Cowie, Songaila, fc Barger 199!: , and 



culation of dust extinction as a function of luminosity using 
the Wang & Heckman scaling, this will underestimate the 
dust correction but not by much (see Fig. ^|). Because we 
have assumed that the slope of the UV luminosity function 
is the same for all of these samples, this is equivalen t to 
adding 0.33 in the log to the values tabulated in Table A2. 
To obtain a "maximal" dust correction, we apply the factor 
of 4.7 to all galaxies in the sample. As we discussed in item 
4 above, this is likely to be an overestimate. Note too that if 
we integrate the new Madau diagram with the fiducial dust 
correction, the total mass of stars pro duced by z = is in 
good agreeme nt with the estimates of Fukugita, Hogan, & 



Peebles (1998 ). However, if we integrate the Madau diagram 
with the maximal dust correction, the total mass in stars is 
a factor of five too large (see Fig. |l^). This suggests that 
something in between the minimal and maximal dust cor- 
rection is probably reasonable (though one could get around 
this conclusion by varying the the IMF from the assumed 
Salpeter shape). 

We plot the results for our fiducial cosmology in Fig. ^. 
For comparison we also plot the results at low redshift from 
recent H-q surveys, corrected for in completeness a n d dust 
extinction by th e original authors ( Gronwall 1995 ; Tresse 
& Maddox 1998). These surveys are quite deep, so incom- 
pleteness corrections are less important, and because Ha 
is emitted at ~6500A, the effects of dust are also less se- 



vere. These results agree well with the Treyer et al. (1998) 
point. There is also a pleasing consistency between the dust- 
corrected, UV-based results at intermediate redshift a nd the 
Far-IR results based on ISO observations at z ~ 0.7 (Flores 



et al. 1999), and finally the high redshift UV-based results 
and the sub-mm results from SCUBA observations at z ~ 3 



( Hughes et al. 199. 



A very different star formation history emerges from 
this "new" Madau plot. Instead of the steep rise from z = 
to z ~ 1.5, there is a more gradual rise, and instead of 
the peak at z ~ 2 and fall-off at higher redshift, there is a 
plateau. 
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